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We discuss possible patterns of electron fractionalization in strongly interacting electron systems. 
A popular possibility is one in which the charge of the electron has been liberated from its Fermi 
statistics. Such a fractionalized phase contains in it the seed of superconductivity. Another pos- 
sibility occurs when the spin of the electron, rather than its charge, is liberated from its Fermi 
statistics. Such a phase contains in it the seed of magnetism, rather than superconductivity. We 
consider models in which both of these phases occur and study possible phase transitions between 
them. We describe other fractionalized phases, distinct from these, in which fractions of the electron 
themselves fractionalize, and discuss the topological characterization of such phases. These ideas are 
illustrated with specific models of p-wave superconductors, Kondo lattices, and coexistence between 
d-wave superconductivity and antiferromagnetism. 
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I. INTRODUCTION 

Electron fractionalization in strongly interacting elec- 
tron systems in dimensions larger than one has been an 
important subject of study since spin-charge separation 
was suggested as a mechanism of high T c superconductiv- 
ity in the cuprates. In particular, it was suggested 
that the electron is splintered into a spin-carrying neutral 
excitation ('spinon') and a charge-carrying spinless exci- 
tation ('holons' or 'chargons'). There have been different 
proposals in regard to this possibility, but the existence 
of such phases in the cuprates is still controversial. 

On the other hand, there exist clear experimental ex- 
amples of phases in the quantum Hall regime of two di- 
mensional electron systems where quantum number frac- 
tionalization has been well established. The low en- 
ergy excitations (quasiparticles) in these two-dimensional 
strongly interacting electron systems carry fractions of 
the quantum numbers of the original electrons. Differ- 
ent quantum Hall liquid states can be characterized by 
different varieties of topological order. The transitions be- 
tween different quantum Hall states can be understood as 
topological-order-changing transitions which occur even 
in the absence of conventional broken symmetries. The 
Hall conductance is but one of the topological quantum 
numbers which characterize a given phase. Another im- 
portant property of a topologically ordered state is the 
ground state degeneracy of the system on higher genus 
manifolds such as tori. For each topologically ordered 



state, there are corresponding sets of characteristic exci- 
tations with different quantum numbers. 

It has become clear |3|Q] that the notion of topologi- 
cal order also provides a precise characterization of spin- 
charge separated and other fractionalized phases in spa- 
tial dimensions higher than one even in situations of zero 
or weak magnetic fields. One of the remarkable features 
of the quantum Hall effect is the enormously rich number 
of exotic phases which display different patterns of frac- 
tionalization of the electron and associated topological 
orders. In view of the similarity between the theoretical 
characterization of quantum Hall states and fractional- 
ized states in zero magnetic field, it is tempting to inves- 
tigate a similar possibility of a variety of fractionalization 
patterns in other strongly correlated systems. We explore 
this possibility in this paper. We describe theoretically a 
few of the several different possible fractionalized phases 
that may exist in various different models of strongly in- 
teracting electron systems. 

Following the introduction of the Schwinger boson de- 
scription of the Heisenberg model of quantum antiferro- 
magnets || , slave fermion j^] formulations of doped anti- 
ferromagnets were introduced. In these formulations, it 
is assumed that the electron decays into a bosonic, spin- 
1/2 spinon and a fermionic, charge-e holon. We will call 
this phase CfSb (Charged /ermion, Spinful &oson). 

On the other hand, a phase with bosonic holons and 
fermionic spinons - which we will call CbSf (Charged 
&oson, Spinful /ermion) - naturally leads to superconduc- 
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tivity through the Bose condensation of bosonic holons in 
the presence of fermionic spinon pairing. Consequently, 
much attention has been focused on the description of 
such a fractionalized phase, especially in the context of 
the slave boson description of the t — J model. The pair- 
ing symmetry of the resulting superconductor is dictated 
by the underlying symmetry of the spinon pairing. 

A Zi gauge theory of fermionic spinons and bosonic 
holons was developed in the context of superconductivity 
in the cuprates (see also [0,0, and Q ) . Spinons and 
holons are coupled by an Ising gauge field. The decon- 
fined phase of this theory corresponds to the CbSf phase. 
Most importantly, the deconfinement-confinement tran- 
sition of spinons and holons occurs through the conden- 
sation of vortices in the Zi gauge field, or visons. In the 
deconfincd phase, the visons exist as gapped excitations; 
when visons condense, the spinons and holons are con- 
fined within electrons. The existence of gapped visons 
is crucial for the robustness of the topological order of 
the deconfined fractionalized phase jLl],|j . Although this 
formalism was introduced in the context of cuprate su- 
perconductivity, it is sufficiently flexible to permit a de- 
scription of other types of fractionalized phases including 
CfSb. 

These ideas have a physical manifestation in the con- 
text of quantum disordered magnets and superconduc- 
tors. In this picture, one visualizes fractionalized states 
in terms of nearby ordered states. In a broken (con- 
tinuous) symmetry state, Goldstone modes can screen 
the associated quantum number(s) Thus it is pos- 

sible for quasiparticles to be stripped of some of their 
quantum numbers. One might imagine that the de- 
struction of order by quantum fluctuations can pre- 
serve this screening of quasiparticle quantum numbers. 
This occurs when those topological defects of the or- 
dered state which braid non-trivially with the quasipar- 
ticles persist as gapped excitations even after the demise 
of the order |13|| . Indeed, this is precisely what hap- 
pens when the state is topologically-ordered. The neu- 
tral, spin-1/2 fcrmion of the CbSf state is viewed as 
the descendent of the Bogoliubov-de Gennes quasipar- 
ticle; the vison, of the hc/2e vortex. When consid- 
ered in the context of spin-triplet superconductors and 
their rich order-parameter structure, this immediately 
suggests exotic phases such as CbSf, CfSb, and even 
a third phase CbSbNf (Charged ftoson, 6pinful ftoson, 
iVeutral /ermion) , in which the charge- and spin-carrying 
excitations are bosons and there is a neutral, spinlcss 
fermionic excitation. Since these superconductors can 
break both charge and spin symmetries - as do states in 
which singlet superconductivity and magnetism coexist - 
one can envision the screening of both quantum numbers 
of a quasiparticle. If (the minimal) topological defects 
in the charge sector survive into a disordered state, then 
this disordered state has neutral, spin-1/2 fermionic ex- 
citations (CbSf ); if topological defects in the spin sec- 



tor survive into a disordered state, then this disordered 
state has charge-e spinless fermionic excitations (CfSb); 
if topological defects in the both sectors survive into a 
disordered state, then this disordered state has neutral, 
spinlcss fermionic excitations (CbSbN /). 

The analysis of quantum dimcr models || and 
resonating-valence-bond Jl|, pT| , ^5[ g round states led to 
conflicting claims that the CbSf^ or CfSb 0,|l8| sce- 
nario is realized in these models. These models have a 
Z2 vortex excitation |T(^,|l7j - which are precisely the vi- 
sons described above - which are relative semions with 
spinons and holons. Thus, a spinon or holon can change 
between bosonic and fermionic statistics by forming a 
bound state with a vison. This begs the question whether 
the CbSf phase discussed in the context of superconduc- 
tivity is the same as the CfSb phase considered in rela- 
tion to magnetism. We reconsider this question in the 
context of recent progress in the understanding of frac- 
tionalized phases described above. One might worry that 
the apparent differences between these phases is an arti- 
fact of the formalisms employed. One might also wonder 
if there are any further fractionalized phases. In this pa- 
per, we discuss the questions raised above using three dif- 
ferent models: p-wave superconductors, Kondo lattices, 
and XY magnets coupled to d-wave superconductors. 

The main results can be summarized as follows. 

1. Both CbSf and CfSb phases can arise in a variety of 
different models. 

2. Upon accepting the possibility of electron fractional- 
ization, one is led to consider a wide variety of fraction- 
alized phases. In the higher-level fractionalized phases, 
electrons can be fractionalized in many different ways. 
For example, spinons and holons can be further frac- 
tionalized. Apart from the CbSf and CfSb phases, we 
discuss two others. One is the CSbNf ( Charge- and .Spin- 
carrying 6oson, ./Veutral fermion) phase, in which the 
electron breaks up into a boson which carries both the 
spin and charge quantum numbers and a neutral fermion. 
This phase is at the first level of fractionalization along 
with the CbSf and CfSb phases. The other is the Cb- 
SbNf phase, in which there exist spin-carrying neutral 
bosons, charge-carrying spinless bosons, and "statistics- 
carrying" neutral spinless fermions. The CbSbNf phase 
is at the second level of fractionalization. In principle, 
higher-level fractionalized phases exist. 

3. We demonstrate the existence of some of these exotic 
phases in the context of the three different systems men- 
tioned above - Kondo lattices, p-wave superconductors, 
and models with both strong spin and d-wave pairing 
fluctuations. For the p-wave superconductor, the four 
fractionalized phases discussed here arise naturally or- 
der parameter has a rich spectrum of topological defects 
which can condense in a variety of ways, thereby giv- 
ing rise to an array of fractionalized non-superconducting 
phases. 
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4. The question of whether CbSf and CfSb are smoothly 
connected to one another or whether they are necessarily 
separated by a phase transition is a subtle and delicate 
issue for reasons that will be discussed at length later. 
While we do not provide a definitive conclusion, we out- 
line a possible scenario in which the distinction between 
CbSf and CfSb is similar to that between liquid and 
gas phases. These phases are separated by a first order 
transition line which terminates at a critical point. In 
principle, one can go around the critical point from one 
phase to the other without encountering a phase transi- 
tion. This scenario is supported by a number of sugges- 
tive (though certainly not conclusive) arguments. 

On the other hand, the transition between the two 
phases can occur through another fractionalized phase 
with a higher-level fractionalization pattern. In this case, 
each transition in the process could be a continuous tran- 
sition. We demonstrate that the transition between CbSf 
phase and CfSb phase can occur through the CbSbNf 
phase. 

5. In order to examine whether one can go from CbSf 
to CfSb through further fractionalized phases like CbS- 
bNf, one can design a gedanken flux trapping experiment 
similar to the one proposed in |l9| . This gedanken ex- 
periment clearly demonstrates the existence of a phase 
boundary between CbSf and CfSb when these phases are 
close to CbSbNf. 

Topological order is robust against local perturbations 
such as impurities. Thus we will concentrate on general 
universal properties of the fractionalized phases. One 
of our goals will be to give a precise characterization of 
these phases which is independent of the underlying mi- 
croscopic models where they may occur. We believe that 
these exotic phases could play a role in the physics of 
3 Hc [^(| and the ruthenates pl|] as well as the cuprates 
p2[ , organic superconductors |23|, heavy- fermion super- 
conductors [^4), spinor Bose-Einstein condensates [ B5J , 
and the crusts of neutron stars p6[ . 

The rest of the paper is organized as follows. In section 
II, we consider a Kondo lattice model and how the CfSb 
fractionalized phase can occur in this model using the 
language of a Z2 gauge theory. Some details are given 
in Appendix A. In section III, we suggest how this anal- 
ysis can be generalized and discuss a hierarchy of frac- 
tionalized phases. Here we provide an overview of our 
results. In section IV, we discuss how this hierarchy can 
be realized in p-wave superconducting systems when the 
superconducting and spin order are quantum disordered. 
This is done using the vortex condensation formalism. 
In Appendix B, the same ideas are shown to apply to an 
XY magnet which is coupled to a d-wave superconductor. 
In section V, the fractionalized phases of section IV are 
further discussed in the framework of a Z2 x Z2 gauge 
theory. In section VI, we consider the question of the 
distinction in principle between the putatively-different 



fractionalized phases constructed in this paper. In Ap- 
pendix C, we give some technical details of an argument 
using Z2 x Z2 gauge theory which supports our picture 
of the phase diagram. In section VII, we show how flux- 
trapping experiments (of the variety suggested by Senthil 
and Fisher (jjj) can be used to shed further light on the 
phase boundaries between these phases and could be used 
to detect them. We conclude in section VII. Appendix D 
contains an aside in which we discuss various novel prop- 
erties of unfractionalizcd phases occurring in the models 
considered in this paper. 

For other perspectives on fractionalization, see 




II. FRACTIONALIZATION IN SPIN MODELS: 
SPIN-STATISTICS SEPARATION 

In principle, there are several possible ways in which 
the electron can fractionalize in a strongly correlated sys- 
tem. In the context of the cuprates, attention has focused 
on the situation in which the electron splinters into two 
separate excitations -a charged spinless boson, and a neu- 
tral spinful fermion. In this case, the charge of the elec- 
tron is liberated from its Fermi statistics. 

In this section, we will briefly discuss another possible 
fractionalization pattern in which the spin, rather than 
the charge, of the electron is liberated from its Fermi 
statistics. The electron splinters into a charged spinless 
fermion, and a spinful boson. As we will see, this phe- 
nomenon also requires the presence of a gapped topo- 
logical Z2 vortex excitation. The issue of whether such 
a fractionalized phase is distinct from one in which the 
charge is liberated from the Fermi statistics is a delicate 
one, and shall be discussed in the Section [vi|. 

To motivate the discussion, consider a "Kondo-lattice" 
model with the Hamiltonian 

H = H t +H K + H ex (1) 
H t = - 2_j trr ' ( c ™<Va + h.c) (2) 

<rr / > 

h k = JkJ2 (^ + 4i c n + h -c) (3) 

r 

H - = J2~^ ( S r S r> +h.c)+ J z S Z r S Z r , (4) 
rr' 

Here the Ci a represent "conduction" electrons with spin a 
at site i. The operators Si are spin operators representing 
magnetic moments localized at the lattice sites. The first 
term is the usual conduction electron hopping, described 
in a tight-binding approximation. The second term is a 
"Kondo" coupling between the conduction electrons and 
the local moments. The third term is an explicit ex- 
change interaction between the local moments. For sim- 
plicity, we have assumed that system only has a U(l) spin 



3 



symmetry for rotations about the z-axis of spin (we will 
comment on situations with full SU(2) spin symmetry 
later). We are interested not so much in establishing the 
exact phase diagram of this particular model - rather our 
main interest here is in establishing the possible existence 
and stability in models of this kind of quantum phases 
where the electron is fractionalized. To that end, we will 
think more generally about a class of models which may 
be obtained from the model above by adding other local 
interactions which share its symmetries. If the system 
is in a quantum phase in which both the symmetry of 
rotations about the ^-direction of spin and the charge 
conservation symmetry is unbroken, the excitations may 
be labelled by their S z and charge (Q) quantum num- 
bers. Clearly, we can visualize two qualitatively differ- 
ent possibilities. First, the system may be in a phase in 
which the excitations are electrons (Q — l,S z — I) or 
composite objects made from electrons (such as, for in- 
stance, a magnon which has Q = 0,S Z = 1). This is a 
conventional phase of the kind familiar from textbooks 
(for instance, a Fermi liquid or a band insulator). On the 
other hand, one could also imagine phases in which there 
are excitations which carry quantum numbers which are 
fractions of those of an electron. The simplest possibility 
(the one we will focus on) is that there are excitations 
which carry S z — 1/2, Q = ("spinons") and others 
which carry S z = 0, Q = 1 ("holons"). In such a phase, 
the electron has been fractionalized. In what follows, we 
will discuss several ways of thinking about such phases. 
Our focus will be on general universal properties of such 
phases. In particular, we will be interested in obtaining 
robust precise characterizations of fractionalized phases 
that are independent of the particular microscopic mod- 
els in which they possibly occur. 

It is extremely instructive to begin by just considering 
the physics of the local moments alone as described by 
the exchange part of the Hamiltonian H ex . This Hamil- 
tonian is clearly invariant under a global spin rotation 
about the z-axis of spin. For technical simplicity, we will 
assume J, J z > 0. The physics of this particular Hamil- 
tonian is well-understood: when J z /J is small, there is 
long range order in S + . When J z /J is large, the system 
breaks translational symmetry with < S z > being larger 
in one sublattice of the square lattice than the other, but 
the U(l) spin rotation symmetry is unbroken. The point 
J z — J can be mapped to the nearest-neighbour antifer- 
romagnetic Heisenberg model with full SU (2) spin sym- 
metry on a bipartite lattice by rotating the spins on one 
sublattice by tt about the z-axis. In the specific case of 
a square lattice (which we assume through out our dis- 
cussion), this is known to develop Neel long range order 
in two spatial dimensions. Our interest here is not so 
much in the properties of this particular Hamiltonian as 
in the properties of an entire class of systems with the 
same symmetry, and with short-ranged interactions be- 
tween the spins. In particular, we will be interested in 



fractionalized phases in which the excitations are spinons 
with quantum number S z = 1/2. To that end, we will re- 
formulate the Hamiltonian directly in terms of "spinon" 
fields which carry spin S z = 1/2. This naturally intro- 
duces a Z 2 gauge symmetry. The result is a theory of 
bosonic spinon fields coupled to a Z 2 gauge field which 
can then be used to analyze the possibility of fractional- 
ized phases and their universal properties. 

We may think of S + , S~ as the creation and destruc- 
tion operators respectively of a hard core boson on the 
sites of the lattice. Specifically, write 5+ = b' sr , S~, = b sr , 
and S z = 1/2 — b\ r b sr . Note that there is half a boson for 
each site on average. Now imagine relaxing the hard-core 
constraint on the bosons, and instead add a term 

|]T(2n r -l) 2 (5) 

r 

at each lattice site. Here n r is the boson number at each 
site. In the limit U —> oo, we recover the spin model 
exactly. For large but finite U however, relaxing the hard- 
core constraint is expected to be innocuous. It is now 
convenient to go to a number-phase representation for 
the bosons: we write b sr ~ e %<Pr with [ip r , n r i] = i5 rr >. For 
simplicity, we also specialize to the limit where J z = 0. 
The Hamiltonian then becomes 

H= -J cos((/J r — <p r >) + ^ ^(1 — 2n r ) 2 (6) 

<rr'> r 

This is clearly closely related to the original spin Hamil- 
tonian in Eqj4[ Now consider a formal change of variables 
which involves splitting the boson operator b sr into two 
pieces: 

b sr = e ivv = z 2 r (7) 
z r = e^= s r e i£ f- (s r = ±1) (8) 

We will refer to the spinon destruction operator. 

Note that with these definitions, both (p r and cj> r are de- 
fined in the interval [0, 2tt). It is also convenient to de- 
fine a number operator for the spinons N r = 2n r which 
is conjugate to <f> r . In terms of the spinon operator, the 
Hamiltonian becomes 

H= ^ -J cos(2</> r - 2<j) r >) + — ^2(N r - l) 2 . (9) 

<rr'> r 

The change of variables above must be supplemented 
with a constraint - clearly the physical Hilbert space con- 
sists only of states where N r is even. Therefore, we need 
to impose the operator constraint (— l)^ = 1 at each 
site of the lattice. Formally this may be implemented 
through the projection operator 

V = Y[V r (10) 

r 

TV = i(i + (-!)*•) (11) 
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Note that [P, H] = 0. It is now convenient to pass to 
a functional integral formulation. We follow Ref. Jf|j27| 
closely to obtain for the partition function 



z 


= E / v ^~ s 


(12) 


s 


= S T + S r + Sb 


(13) 




= E JrCTrr COs(<f» r ,T+e - (j> rT ) 
r,r 


(14) 


Sf 


= e 2_, J cos(20 rr — 2<ft r > T ) 


(15) 




<rr'>r 




rr 


= ±1 may be interpreted as 


the time- 



component of a Z 2 gauge field that imposes the constraint 
on the Hilbert space, and e is the lattice spacing along 
the time direction. The constant J T is determined by the 
original interaction strength U. The term in the action 
S r involving the spatial coupling may be decoupled by a 
Hubbard-Stratanovich transformation: 

£ -S r = J Dxe -^Y, <Trl> , T X rr> {r) 2 +2eJ Xrrl (r){zl(r)z Tl (r)+c.c) 

(16) 

Here Xrr'( T ) i s a real- valued field. We have omitted an 
unimportant overall constant. 

We now proceed exactly as in Ref. and replace 

the integral over the continuous variable x by a sum over 
a discrete field a rr >{r) = ±1. As discussed in Ref. RCT, 
this approximation respects all the symmetries of the ac- 
tion, and is expected to be innocuous. The resulting 
partition function becomes 

z = E / v ^' s 

S = S S + S B (18) 

S's = - E Ji i ai 3 C0S (<^ ~ <Aj) ( 19 ) 

<ij> 

Here the i,j label the sites of a space-time lattice in three 
dimensions. The constants Jjj = J T for temporal links, 
and equals eJ for spatial links. Sb is the Berry phase 
action 

Sb = y E t 1 -^) ( 2 °) 

i,j=i-\-T 

Note that the action (18)-(20) respects all the symmetries 
of the original model. The discrete field = ±1 may 
be viewed as a Z 2 gauge field. At this stage, this field 
does not have any dynamics. However, it is natural to 
expect that upon coarse-graining, some dynamics will be 
generated. The simplest such term allowed by symmetry 
is 

S K = -KJ2U^- ( 21 ) 

□ □ 



We will therefore consider the full action 

S = S s + S K + S B (22) 

What we have achieved so far is an approximate re- 
formulation of spin models with XXZ symmetry. This 
reformulation is extremely useful to explore the various 
possible allowed phases in such models. However, the 
approximations made in obtaining this reformulation are 
severe enough that it is not easy to see which one of 
these allowed phases will be obtained in any particular 
microscopic model. 

Consider the possible phases when the parameter K 
is very large. When K — oo, the Z 2 flux through each 
plaquette is constrained to be one. We may then choose 
a gauge in which Oij = 1 on every link. In this limit 
therefore, the action reduces to 

S = - E J ij COS (0i - 4>j) (23) 
<ij> 

This simply describes a quantum XY model in two spa- 
tial dimensions. Note that the Berry phase term simply 
vanishes when all the Oij = 1. There clearly are two pos- 
sible phases - an XY ordered phase in which = e 1 ^ 1 
has condensed, and a paramagnetic phase in which the 
excitations created by Zi are gapped. Note that these ex- 
citations in the paramagnetic phase carry spin S z = 1/2. 
Thus, the spin has been fractionalizcd in this phase. 

Now consider moving away from the limit K = oo 
by making K large but finite. For finite K as can be 
seen from the arguments advanced in Ref. fit], the XY 
ordered phase where the spinon field has condensed is 
indistinguishable from a conventional XY ordered XXZ 
magnet. The paramagnetic phase in which the spinons 
are uncondensed and deconfined survives for large but fi- 
nite K . When K is finite, it becomes clear that this phase 
has another distinct excitation which carries the flux of 
the Zi gauge field. This Z^ vortex - dubbed the vison - 
does not carry any physical spin, and has an energy gap 
of order K for large K. It has the important property 
that when a spinon is taken around it, the wavefunction 
of the system acquires a phase of ir. 

Upon decreasing K, at some critical value, the vi- 
son gap goes to zero. For smaller K, the visons con- 
dense leading to confinement of the spinons. The result- 
ing phase is a conventional quantum paramagnet with 
gapped S z = 1 excitations. In this phase, the Berry 
phase term becomes important and leads to a breaking 
of translational symmetry - the paramagnet is therefore 
expected to develop spin-Peierls order. We will not dis- 
cuss such conventional phases very much in this paper. 

Much further insight into the physics of the fraction- 
alized phases may be obtained by the following consider- 
ations. We begin by first considering ordered phases in 
which the symmetry of rotations about the z direction of 
spin has been broken spontaneously. For simplicity, we 
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consider a phase in which the spins have all lined up along 
some direction in the xy plane. The general properties 
of such a phase are well-known. There are two distinct 
kinds of excitations. First, there is a gapless spin- wave 
mode with linear dispersion. Apart from these, there are 
also topological vortex excitations. On moving along any 
circuit that encloses a vortex, the direction of the spin in 
the xy plane winds by an integer multiple of 2-7T. This 
integer winding number - the vorticity - is conserved, 
and may be used to label the spectrum of excited states. 
States with different total vorticity belong to different 
topological sectors and are not mixed by the dynamics 
generated by the Hamiltonian. Note that in this ordered 
phase we can no longer label states by their S z quantum 
number. 

These familiar properties of the XY ordered phase 
must be contrasted with those of the quantum param- 
agnet. First consider a conventional paramagnet (i.e. 
one with no fractionalization) . Clearly in this phase S z 
is conserved, and is a good quantum number to label the 
excitation spectrum. On the other hand, the vorticity 
loses its meaning in the paramagnetic phase, and is no 
longer a good quantum number. This suggests that one 
may view the paramagnet as a phase in which the vortex 
excitations have themselves condensed. Condensation of 
the vortices implies that the vorticity is no longer a good 
quantum number (just like condensation of spin implies 
that S z is no longer a good quantum number). Indeed, 
these observations may be formalized precisely by means 
of a duality transformation which reformulates the sys- 
tem in terms of the vortex fields rather than the spins. In 
this dual formulation, the paramagnet is described as a 
vortex condensate, and the XY-ordered phase as a vor- 
tex insulator (in which the vortices are gapped). The 
physical excitations of the paramagnet which carry the 
S z quantum number appear as dual flux tubes of the 
vortex condensate in this language. 

How are we to view the fractionalized quantum param- 
agnet in this dual language? As the phase in question is 
a paramagnet, it is clear that the vorticity has no mean- 
ing, implying that the vortices must have condensed. As 
pointed out in Ref. JD|, we may view the fractionalized 
phase as a condensate of paired vortices. This has the 
immediate consequence of halving the dual flux tube, i.e 
of fractionalizing S z as required. Furthermore, note that 
the unpaired (and uncondensed) single vortex is still an 
excitation in the system. Its vorticity is screened by the 
(double strength) vortex condensate as is required in the 
paramagnet. However, its parity is still a good quantum 
number. Thus the unpaired vortex, though a legitimate 
excitation of the fractionalized paramagnet, carries only 
a Z2 quantum number - it is clear that it is the vison 
excitation discussed previously. 

The discussion above provides a description of a frac- 
tionalized quantum paramagnet in the context of spin 
models with XXZ symmetry. We now return to the full 



model which includes coupling to the "conduction" elec- 
tron degrees of freedom. As above, we first replace the 
operator S~ in the Kondo coupling at each site by the 



boson operators b s 



(and similarly for S^). The 



Kondo coupling term then becomes 



H K = 4^ (&J r 4jC rT + h.c. 

r 

= 4^ [z^cl^ + h.c 



(24) 
(25) 



In going to the second equation, we have introduced the 
spinon operators z r defined in (Q). The Kondo coupling 
can be further simplified by another change of variables 



T}r1 — Z r C r ^ 
Vrl = 4 c rl 



(26) 
(27) 



We will call the ^-operators the holon operators. In terms 
of the holons, the Kondo coupling becomes 



H, 



K 



h.c 



(28) 



Note that the holons are actually spinless charge e fields 
despite the presence of the label t,j- This is obvious 
from their definition in terms of the spinon and electron 
operators above: the holon operators do not transform 
under spin rotations about the z-axis. Explicitly, the 
Kondo term mixes up and down holons so that their label 
I, J, is changed by the dynamics. Therefore their spin 
label has no great significance, and they are correctly 
viewed as spinless fermions. We may use the following 
physical picture: the Kondo spins screen the spin of the 
conduction electrons. 

Under these changes of variables, the electron hopping 
term becomes 



H, = - 



E 

<rr'> 



h.c 



(29) 



We now make approximations very similar to those used 
above for the exchange part of the Hamiltonian. They al- 
low us to reformulate the system in terms of the spinons, 
holons and a Z2 gauge field. Some of the details are 
oulincd in the Appendix. The resulting action can es- 
sentially be guessed on symmetry grounds, and takes the 
form: 



S = S C + S S + S B + S k (30) 
Sc = - V VijtijiViiVjl + vliVjl + cc.) (31) 



E 

<ij> 



The other terms of the action are as given before. 



G 



Following the discussion above, for large K, we expect 
to have a phase in which the holons and spinons are lib- 
erated from each other. In such a phase, the electron is 
fractionalized. However, in contrast to the fractionalized 
phase that is most popular in the context of the cuprates, 
here the spin of the electron has been liberated from its 
Fermi statistics. Are these two phases actually the same? 
We will address this issue in subsequent sections. 

Though we have based our discussion on models with 
XXZ symmetry, we expect the fractionalized quantum 
paramagnetic phase to exist even in systems with full 
SU(2) spin symmetry. Indeed, in the context of frus- 
trated quantum Sp(ti) spin models in the large-n limit, 
Read and Sachdev || have argued for the stability of 
fractionalized paramagnetic phases with properties simi- 
lar to that discussed above. 

III. A HIERARCHY OF FRACTIONALIZED 
PHASES 

In Section [n|, we primarily discussed fractionalized 
phases in which the electron splinters into a spin-1/2 
neutral boson and a charged spinless fermion. For future 
convenience, we will refer to this as the CfSb (charged 
fermion, spinful boson) phase. This is to be contrasted 
with the fractionalized phases which are popular in the 
context of cuprate physics in which the electron splinters 
into a spin-1/2 neutral fermion and a charged spinless 
boson (see also Sections IV and V). We will refer to this 
as the CbSf phase. In both cases, there is, in addition, a 
Z2 vortex excitation (the vison) such that taking either 
the holon or spinon around it produces a phase change 

of 71". 

Having accepted the possibility of quantum number 
fractionalization, one can imagine a wide variety of pos- 
sible phases apart from the two mentioned above. In 
particular, one may consider exotic possibilities where 
the fractions of the electron in any given fractionalized 
phase themselves fractionalize. Such phases may be con- 
sidered to have a higher level of fractionalization. To see 
how these may be described in the same kind of formu- 
lation as discussed in this section, consider the following 
action: 

S = S f -\~ S c ~\- S s ~\~ Sen- 

Sf = z2 a ^'j ' ''../•■■/.•" +t A ^(Mfji - (T- > l)) + c 

ij,a 

^ 1 fia fia 1 

ia 

S C = ~Y^ tijTij (b* ci b cj + c.c.) , 
y 

Ss = - ^t^CTij^Zj +c.c.) , 
ij 

□ □ □ 



Here b c is a charge e spinless boson and z is a spin-1/2 
chargeless boson. The / field represents a spinless, neu- 
tral fermion (the spin index is just a label with no spe- 
cial significance). The ciy and Ty are two independent 
Zi gauge fields. The physical electron a a = b C iZifi a . 
Clearly if the r field is confining, the f a and b c get con- 
fined to form a fermionic holon - we then recover the ac- 
tion discussed earlier in this section. On the other hand, 
if the a field is confining, the Fermi statistics gets glued 
to the spinon (zj) - the resulting theory is essentially that 
introduced in Ref. B in the context of cuprate physics 
and involves bosonic holons and fermionic spinons cou- 
pled to a Z2 gauge field. If both gauge fields a and r are 
deconfining however, we have an exotic phase in which 
the fields b c , z, f are all liberated. This phase will also 
have two distinct vison excitations corresponding to the 
fluxes of the two Z 2 gauge fields. We may view this phase 
as a higher-level fractionalized phase as compared to the 
one discussed in Ref. or that discussed earlier in this 
section. The connection between various fractionalized 
phases is shown in Figure 0. We use symbols b a and f a 
to label bosons and fermions that carry quantum num- 
bers a = n, c, s, cs (neutral, charge, spin, charge and spin) 
and show the existence of appropriate Zi vortices in each 
phase (for more details see Section V). 

In the sections which follow, we will show how an ef- 
fective action such as that of Eq. [32] can arise in the 
context of p-wave superconducting systems and systems 
which feature interplay between magnetism and super- 
conductivity. 
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IV. FRACTIONALIZATION OF ELECTRON 
QUANTUM NUMBERS WITH P-WAVE 
PAIRING 

A. Order Parameters and Symmetries 

Spin-triplet superconductors and their rich order- 
parameter structure offer the prospect of various exotic 
phases. Since they break both charge and spin sym- 
metries, triplet superconductors exhibit features of both 
singlet superconductors and of spin models. In partic- 
ular, we can envision the restoration of the the £/c(l) 
charge symmetry by quantum fluctuations, thereby re- 
sulting in a spin-triplet insulating state. Alternatively, 
the spin symmetry (we will make the simplifying assump- 
tion that the system has only an easy-plane Us {I) spin 
symmetry) can be restored, resulting in a spin-singlet 
superconducting state. Finally, both symmetries can be 
restored, leading to a singlet insulating state. We believe 
that the gapped, symmetry-restored states will not be 
very sensitive to the precise symmetry of the spin sec- 
tor, so we believe that our results apply to systems with 
full SU (2) spin symmetry as well. In particular, when 
the symmetry is increased (while keeping the size of the 
representation fixed), fluctuations are enhanced, and a 
system is more likely to be in a disordered state. In or- 
der for these symmetries to be restored separately, it will 
be necessary, as we discuss below, for a type of topo- 
logical ordering to occur. This topological ordering is 
essentially spin-charge separation of the charge 2e, spin- 
triplet Cooper pairs. Depending on the way in which the 
symmetries are restored, it is possible for further topo- 
logical ordering to take place, in which case the quan- 
tum disordered states may support excitations with ex- 
otic quantum numbers. In such states, the spin and/or 
charge of the quasiparticles is screened by the Goldstone 
modes (which are themselves separated from each other 
by the higher-level topological ordering) . As we describe 
in this paper, there are no fewer than nine phase which 
can result in this way. 

To be concrete, let us consider the following p-wave 
superconducting state of electrons on a square lattice: 

A a p = A e %v (cosOa^p + isin9 5 a p) sin k y a (33) 

This is the most general unitary triplet state in 2D |2C| l 
if we assume that there is only the U(l) spin symmetry 
of rotations about the z-axis, rather than the full SU(2). 
In (|33l), only A^ and A^ are non-zero. The lower sym- 
metry could be the result of spin-orbit coupling. The 
symmetry-breaking pattern associated with this order 
parameter is: £/c(l) x Us(l) x Z? 4 — > Z 2 x Z 2 X D 2 . The 
t/c(l) charge symmetry is broken to Z 2 by the conden- 
sation of a charge 2e order parameter. The Us(l) spin- 
rotational symmetry is completely broken. The square 
lattice point group, D4, is broken to D 2 by the orbital 



symmetry of A. Finally, there is an additional Z 2 since 
the order parameter is left invariant by if — » ip + n, 
6 — > 6 + 7r. As we discuss later, this can be understood 
as a Z 2 gauge symmetry. From e lip and e* e we can con- 
struct the following ^-invariant order parameters whose 
presence or absence characterizes the phases which we 
consider. In the absence of the triplet p-wave supercon- 
ducting order parameter ( |33| ) , we can characterize states 
by the charge-4e order parameter, 

and the spin nematic order parameter, 

Q = cos26> (35) 

These order parameters define the following quantum- 
disordered states of triplet p-wave superconductors. 

• Charge-4e singlet superconductor: A 4e ^ 0, Q = 0. 

• Charge-4e nematic superconductor: A 4e ^ 0, Q ^ 
0. 

• Spin-nematic insulator: A 4e = 0, Q ^ 0. 

• Spin-singlet insulator: A 4e = 0, Q = 0. 

B. Topological Defects 

The quantum-disordered and topologically-ordered 
states which we will consider can be understood in terms 
of the condensation or suppression of various topological 
excitations. The most basic and fundamental topological 
excitation is a composite formed of a flux hc/4e vortex 
together with a tt disclination 



composite = — ~~ x 



FIG. 2. 7r-disclination - hc/Ae vortex composite 

Along a circuit about such an excitation, both ip and 9 
wind by tt so that any ^-invariant combination is single- 
valued. If such an excitation is at the origin, and r, (f> are 
polar coordinates in the plane, then the order parameter 
is of the form: 

A a0 (r,(j>) =A(r) e 1 ^' 2 x 

^cos ^ a z af} + ie 2 sin ^ S af ^j sin k y a (36) 

where A(0) = and A(cx)) = Ao- The flux is into or out 
of the plane, respectively, for t\ = ±1; the spins wind 



8 



clockwise or counter-clockwise, respectively, for ti = ±1. 
It is instructive to write this as: 



d 2 x dr x 



A TT (r,</> 
A u (r,$ = -A(r) e 



A(r) e"+ sin£; y a 



sin k y a 



(37) 



<9 r - A c t t z - v F r z id x + v F A% 
- A%o z + v F A a x a z T z 
+ i sin 6 crO 



X (42) 



where e± = (ei ± (-2)j2. Hence, 7r-disclination - hc/4e 
vortex composites are vortices in Aff or Aj^ alone. 

These excitations can be combined to form an hc/2e 
vortex which is non-trivial in the charge sector but trivial 
in the spin sector. 



s = ± and x nas a particle-hole index, acted on by Pauli 
matrices r; and a spin index, acted on by Pauli matrices 
a. 



A a| a(r, 4>) — A(r) (cos^oer^ + isinflo'W) sin fc y a (38) Xaa(fc) 

with constant #o- Alternatively, we can form merons, 
which are trivial in the charge sector but not the spin 
sector. 

A Q/3 (r, 4>) — A(r) e llf0 (cosf/xr^ -I- « sm^Sap) smk y a (39) 

with constant (po. Finally, there are various composites 
formed from the above. A composite formed by n hc/2e 
vortices and m merons takes the form 







c- 

fep+feT 

c} - 

— kp—kl 


X21 




Xl2 




k F +ki 




X22 _ 




-C+ - 

— fcjT — 



(43) 



A a/ 3(r 



A e 



irup 



\(\)<j z a n + ism.m<j>5 a p) s'm k y a (40) 



If flux hc/Ae vortex-7r disclination composites con- 
dense, then t/c(l) an d Us{l) are restored. The sys- 
tem will be in a singlet insulating state and all excita- 
tions will have conventional quantum numbers. If, on 
the other hand, hc/Ae vortex- 7r disclination composites 
are gapped and only complexes consisting of multiples 
of hc/4e vortex- 7r disclinations (e.g. n hc/2e-m meron 
composites) are condensed, then quantum number sepa- 
ration is possible. If complexes consisting of a multiple of 
four hc/4e vortex- tt disclinations condense, then we will 
have the various versions of quantum number separation 
summarized in Figures |5| and ^. 



C. Quantum Number Separation 

The effective action of a p-wave superconductor may 
be written in the form: 




FIG. 3. Order parameter for a sink y p-wave superconduc- 
tor. Gapless excitations exist at fcj? = (±&f,0). 



In action 



2|) we have included the electromagnetic 

(7 
I' 



field A c „ and spin vector potential A a which couple to 



the conserved electric and S z currents. 

When hc/Ae vortex- ir disclination composites are 
gapped, the Zi symmetry yj— »y> + 7r,0— »0 + 7r plays no 



role, and the other terms in (41) may be written in the 
form. 



5*tot — Sf + S c + S a 



(41) 



and 



where Sf is the action for the fermionic quasiparticles 
and their interactions with the Goldstone modes, and S c 
and S a are the actions for the charge and spin Goldstone 
modes. 

Depending on the topology of the Fermi surface, the 
low-energy spectrum of a p-wave superconductor may in- 
clude gapless fermionic quasiparticles. Let us assume 
that the topology is such that the gap has nodes on the 
Fermi surface. Focusing on the nodes, as shown in figure 
0. We linearize the action: 



Sa = \po J d 2 xdr (dJ-A°) 



(44) 



(45) 



The conserved electric and S z currents are given by: 
8S tot 



J x SAT 



(46) 



Conservation of charge and the z-component of spin re- 
quire 



W = 



(47) 
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The interactions between the Goldstone fields and the 
quasiparticles are highly non- linear in (0). This inter- 
action can be made more tractable, following p3|, if we 
define new fermion fields ip: 



_ J?T*/2 p i9a*/2 



1> 



(48) 



With this change of variables, we have defined a neu- 
tral, spinless fermion, ip, which is governed by the action: 

Sf = J d 2 xdr y tp^[d T - v F r z id x - v a t x <j x (id y )]ip 

+ ~ ^[ T z d T <p - 2A c t t z - v F 8 x p + 2v F A%]il> 
+ ~ ip\a z d T 9 - 2A°t z - v F T z a z d x 9 + 2v F T z a z A° E }iP 

The couplings between the Goldstone modes and the 
quasiparticles are now either trilinear or biquadratic 

S f = S° f + 1 j d 2 xdr (jS(d T (p - 2A%) 

+ J°(d x p-2A c x ) + Jg(d T 6-2A°) 

+ r x {d x 0-2Al)\ (50) 



with 



J§ = 1p^T Z 1p 

J ff = ^o z i> 



-V F 1p'tp 



J x = -v F ip^a z T z ip 



(51) 



The price that must be paid is that the change of vari- 
ables (ji^) is not single- valued about a topological defect. 
In particular, the charge part, exp(ipr z /2), is double- 
valued under transport about a flux hc/2e vortex since p 
winds by 2n, while the spin part, exp(i9a z /2), is double- 
valued under transport about a meron since 9 winds by 
2tt. 

As we will see below, "0s are weakly-coupled quasipar- 
ticles in those quantum disordered phases in which flux 
hc/2e vortices and merons are gapped. 



D. Defect Condensation 

Defect condensation is now implemented with dual rep- 
resentations for the order parameters. j2^,|l^,^,Q . In 
the dual description of the XY model [ p4| , the ordering 
field, ip, is replaced by a gauge field which parametrizes 
the total current, together with a vortex field which ac- 
counts for the singularities of p. 

We use the conservation of charge to define the dual 
gauge field 



e^xd v al = 4 ot c = Pc (d^p - AD + j; 



(52) 



with from (|51j) , and introduce the vortex current 



2tt 



£,-Lv\d v d\ip 



(53) 



which is not vanishing for a multivalued ip. With the 
last two equations we can relate the vortex current to 
the dual gauge field and quasiparticle current 



il = e^A<9„ [p c 1 e Xa pd a a% + A c x -p c 1 J£] 



(54) 



Now a dual action for the charged degrees of freedom 
is easily constructed by requiring that its equations of 
motion reproduce (|54|). 



S Dual = SGL{^v,a C ^) + 

J drd 2 x (J-(f^ f + a^dAAj 



Pc 



wj^ere 



(4' 



S G L[$,<g - J drd 2 x (^|(^-m M )$| 2 + V($)) (56) 

and /° = d v a'j 1 — d^a c v . The field may be thought of 
as a vortex creation field. The vortex current is given by 



J I, 



Pd 
2 



(57) 



An identical construction is now used for 9 with p c 
replaced by p a and J\ by J%: 



S Dual = S G l(®t 



drd 2 



^-{f^f + ale, vX d v {Ax-^-Jx) 

^Po Pc 



where 3?}^ is a meron creation operator. Analogous topo- 
logical objects in the spin sector have been discussed in 
the context of quantum Hall systems j3^-|37j and quan- 
tum antiferromagnets |l^ , |3^ ^^0|] . 

Actions ([sit ) and (|5§| ) need to be supplemented by 
Chern-Simons gauge fields which enforce the minus sign 
which is acquired when a ip encircles a flux hc/2e vortex 
or a meron |2^^l[. With these additions, we obtain the 
following dual action 



~>Dual = SGL{^v,a C 



dr d 2 x 



2p, 



if') 2 + aU^ x d„(Al 



-A 



2al l e^ x d„a s x 1 + a 1 ^ 



dr d 2 x 



2 Pa 



{fLY + <^xd v {Al 



Pa 



-Jx) 



+ 2ale^ x d v af + a 2 



where a 1 ^ 2 and a* 1,2 are the gauge fields that perform 
the flux attachement and enforce the minus sign. 

With this action in hand, we can now address the quan- 
tum disordered phases and quantum number separation. 
In essence, there are three quantum numbers: charge, 
spin, and electron number modulo 2. These can separate 
in a variety of patterns. 
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• If (<&„) 7^ 0, flux hc/2e vortices condense. The 
Meissner effect associated with this condensate im- 
poses 

a£ + < = (60) 

Recalling that e a pd a a c p — Jq°* c (a, (3 — x,y) is the 
charge density and tapdaO^ = Jg is the quasipar- 
ticle density, we conclude that in this phase charge 
is attached to the ips. 

• If ($ m ) 7^ 0, morons condense, and the Meissner 
effect associated with this condensate imposes 

<+< = C) (61) 

As e a pd a a'p = Jq°* a is the local spin density and 
t a f3d a a S {3 = Jq we find that spin is attached to the 
-0s. All the fcrmions carry spin. 

• If <$ m > = 0, <$„) - 0, but ($„$ m ) ^ hc/2e 
vortex - meron composites condense. The Meissner 
effect associated with this condensate imposes 

*% + < = (62) 

In other words, spin and charge are confined, but 
the fermion -0 carries neither since ip does not ac- 
quire any phase upon encircling this composite ob- 
ject, as evinced by the fact that & v $ m is not cou- 
pled to statistical gauge fields. 

• The condensation of other composites, such as $^ 
(i.e. skyrmions), <I> 2 , etc. does not cause the con- 
finement of any quantum numbers. 

E. Exotic Phases 

The order parameter classification discussed after 
equation ( |35| ) is incomplete; those states can occur in 
several varieties, classified by the allowed quantum num- 
bers f|. 

Charge-Ae singlet superconductors, A 4e ^ 0, Q = 0. 

IA) If (3> m ) 7^ 0, then the fermionic excitations carry 
spin 1/2. 

IB) However, if ($ m ) = but (($ m ) 2 ) ^ then the 0s 
are spinless. Note that the charge quantum number of 
the fermionic excitation is not really well-defined since 
U(l) is broken in the superconducting state; said differ- 
ently, the fermion can always exchange charge with the 
condensate. 

Spin-triplet insulator, A 4e = 0, Q ^ 0. 



2B) If {$.„} = but ({<$> v f) + 0, then the 0s are neutral. 
As in the previous case, the spin quantum number of the 
■0s is not well-defined. 

Spin-singlet insulator, A 4e = 0, Q = 0. 

3A) If ($ t) ) ^ 0, ($ m ) j= 0, then the 0s carry spin 1/2 
and charge e. Phase CSf. 

3B) If ($ m > = but (($ m ) 2 } ^ while ($.„> ^ 0, then 
the 0s are charge e, spinless fermionic excitations: CfSb 
Phase. 

3C) If ($„) = but ((<I\,) 2 } + while ($ m ) ^ 0, then 
the 0s are neutral, spin 1/2 fermionic excitations: CbSf 
Phase. 

3D) If ($„$ m ) 0, then the 0s are neutral, spin- 
less fermionic excitations, but spin and charge are con- 
fined into a bosonic spin 1/2, charge e excitation: CSbNf 
Phase. 

3E) Finally, if ($„) = but ((<5> v f) ^ and ($ m ) = 
but (('I'm) 2 ) 7^ 0, then the -0s are neutral, spinless 
fermionic excitations. Bosonic charge e excitations, e lip / 2 , 
and bosonic spin 1/2 excitations, e l8 / 2 are also liberated: 
CbSbNf Phase. 

To summarize, we have the following phases with ex- 
otic quantum numbers. 

• A charge-4e singlet superconductor with spinless 
fermionic excitations. 

• A spin-triplet insulator with neutral fermionic ex- 
citations. 

• Spin-singlet insulators with (a) charge e spinless 
fermions and spin 1/2 neutral bosons; (b) spin 1/2 
neutral fermions and spinless charge e bosons; (c) 
neutral spinless fermions, bosonic charge e spinless 
excitations, and bosonic spin 1/2 neutral excita- 
tions; or (d) neutral spinless fermionic excitations 
and bosonic charge e spin 1/2 excitations. 

These result are summarized in the following diagrams 
that describe various phases that can result from quan- 
tum disordering a p-wave superconductor. 



2A) If {$.„) 7^ 0, the 0s carry charge e. 
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Quantum disordered 
p-wave superconductor 



Charge 4e singlet 
superconductor 

A 4e >^0 




Merons condensed 
spinons confined 




Merons gapped 
spinons deconfined 
< cj> > = PTr 



— vortices condensed 

:e 

chargons confined 
<<J>> / Wa 



he 

2e 



vortices gapped 



chargons deconfined 



FIG. 4. 
ductor. 



Phases of quantum-disordered p-wave supercon- 



Spin- singlet 
insulator 

<A 4e >=0 

<Q >=0 



he 



2c 



vortices condensed 



merons gapped 
chargons confined 
spinons deconfined 

<V = o 

<<£><!> > : 



he 

2~ vortices gapped 
merons condensed 
chargons deconfined 
spinons confined 
<4>> = 








vortices and merons condensed 
chargons, spinons confined 
<<t> v >^0 



he 

2~ vortex x meron 
composites condensed 
chargons and spinons bound 
deconfined from fermions 

< 4> v > =0 
<d> ra > = 

< *v* m > ^0 [3D 



r- vortices and merons gapped 
chargons, spinons deconfined 
<*„> =0 
< O > = 

m 



FIG. 5. Phases of quantum-disordered p-wave supercon- 
ductor. Note that the phase in which hc/2e vortices 
and merons are condensed may be described as having 
hc/4e- vortex- n composites condensed. 

Scenario proposed in this section for quantum number 
separation in p-wave superconductors may apply to other 
systems, provided that they acquire non-trivial topologi- 
cal order in the spin and charge sectors, or in the language 
of this section when they have sufficiently strong quan- 
tum fluctuations of spin and charge degrees of freedom 
simultaneously. In Appendix ^ we show that quantum 
disordered <i-wave superconductor with easy-plane anti- 
ferromagnetic fluctuations may be treated in the same 
way as we treated p-wave superconductors in this sec- 
tion. 



V. Z 2 x Z 2 LATTICE GAUGE THEORY. 

In this section, we derive a Zi x Zi gauge theory rep- 
resentation of a model which gives rise to local p-wave 
superconducting fluctuations. In addition to the super- 
conducting state, we find the exotic phases discussed in 
the previous section. These have a simple description 
as the various deconfining phases of the Zi x Zi gauge 
theory. Readers who are uninterested in the technical 
details of our derivation may skip directly to equations 
(105), (|106[), and the subsequent discussion. 



A. General Formalism. 

We consider the following Hamiltonian that describes 
the equal spin pairing state of a p-wave superconductor. 



H = H t + H u +H v + H A 



(63) 



with 



H 



~t C ta C r'a + k.C. , 

rr' ,a 

H u = uJ2(N r - No) 2 , 

r 

H v = v^2(M r f , 

r 



(64) 



where AjJ, and A**, represent the order parameter fields 
for the Cooper pairs with spin up-up and down-down 
pairs respectively. Here a =T, J, is the spin index. The 
term proportional to u represents the on-site Coulomb 
repulsion. N r is the total number operator of electrons 
at the site r, Nq is the average electron number per site. 
M r is the z-component of the total spin operator. At 
equilibrium, 

I ATT I = |AU|. Note that there are two in- 
dependent phases associated with ATT ana ; ^11 _ \ye can 
rewrite H A as 

Ha = A J]] a rr >[e liPr1 c r fC r >i + e liPrl c r ic r 'i] + h.c. , (65) 

rr' 

where A = |ATT| = |ATT| and a rr < is the form factor that 
gives rise to the particular p-wave symmetry. 

The fields tp r ^ and (p r ± are canonically conjugate to the 
Cooper pair number operators of up-up and down-down 
Cooper pairs, n r -\ and n r ^: 



[V>rh n r't] = 



= iS r 



(66) 



The conserved charge densities for the electrons with spin 
t and I are given by 

iV r I = 2n r f -|- p r | 
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N rl = 2n rl + p rl , 



(67) 



(4) 



t\2 _ Jlfr 



(76) 



where p ra = c\ a c ra is the quasiparticle number, which is 
not equal to the electron number. It is useful to remind 
the readers that the Hamiltonian (64) does not conserve 
the quasiparticle number, since it contains terms that 
annihilate a pair of them and create a Cooper pair. Only 
the total number of electrons of a given spin, given by 
equation (67), is conserved. This may be formulated as a 
conservation of the total charge and z-component of the 
total spin. 



N r = iV rT + N r l, M r = N rt - 



(68) 



Let us define boson operators, b ra , which carry charge 
e and spin a =|, j: 



&L = tee**™/ 2 = 



(69) 



where i" = ±1 are Ising variables and <p ra are defined in 
the interval zero to 2ir. Note that the squares of and 

&L create the spin up-up and down-down Cooper pairs 
via the following relation. 



2 = e i<Pr 



(70) 



One can also see that the canonical conjugates of <pra are 
the total densities of electrons with spin f and {. They 
satisfy the following commutation relations. 



(71) 



Similarly, the following commutation relations are also 
satisfied. 



(72) 



where <^> rc = (4> r i + r jJ/2 and rs = (0 r | — (f> r i)/2. At 
this stage, it is useful to define the fermion operators, 

fret 5 as follows. 



(73) 



Note that /t Q creates spinless neutral fermions due to 
the fact that 6j Q carries both the charge and spin of the 
electrons. 

It is also useful to define (p rc and tp rs as follows. 



(74) 



Note that there is a Zi symmetry associated with these 
definitions of phase variables; ip rc — > <y9 rc + 7r and <y9 rs — > 
(p rs + 7r do not change e* v ' T and e Jl/3 ' ^ . Now we define 
boson operators, b\ and z}., as 



= s r e" rrs/ " = e" v ' rs . (75) 



Here t r = ±1 and s r = ±1 are Ising variables. Note that 
these operators satisfy the following identities. 



Note also that &L and can be rewritten as 
frt i,t 2 t /,t i,t z 

Now the total Hamiltonian can be written as 
H = H t +H u +H v + H A 
with 

+ /i.e. 

i?A = A a^/ (blb r >z r z r >f r if r ^ + blb r ,z r zl,f rl f r , L ] 



(77) 



(78) 



/i.e. 



(79) 



The Hamiltonian is invariant under the following local 
transformations: 

i) Z^r • b r > b ri fra * fra 
fi) Z^a • %r * %r] fra * fra: 
iii) Z^a • ^-r ^ i ^ra * ^ra 

Only two of these transformations are independent, any 
one of them can be represented as a product of the 
other two. Together they form Zi x Z2 gauge symme- 
try, that has three Z2 subgroups as reflected in three 
possible transformations above. These subgroups are dis- 
tinct, but not independent. Z2 X Z2 local gauge symme- 
try is a consequence of the redundancy in the enlarged 
Hilbert space of f ra , b r , and z r . There is a further redun- 
dancy in our description in terms of b r and z r because 
b r — > ib r , z r — > — i z r also leaves all physical quanti- 
ties invariant. This identification allows for the existence 
of flux 7r-disclination-/ic/4e vortex composites which we 
discussed in subsection IV B. As before, we assume that 



these topological defects are gapped so that we can safely 
ignore this identification and take <p rc and 4> rs as defined 
from [0,2tt). 

In order to get the correct Hilbert space of the elec- 
trons, we have to impose two constraints at each site. 



N r + p r | + Pri = even number, 
M r + p^ — Pri = even number . 



(80) 



These can be written as 

(_l\N r +p rr +p rl _ ^ l_l\M r +p rl - Prl _ j _ egy 

The constraints can be implemented in the path integral 
representation of the partition function using the follow- 
ing projection operators. 



^Pc — Pre: P s I I Pr 



(82) 
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with 



V rc = -[l + (-l) Ar "+^T+^i] 



T r )(N r +p rt +p r i) 



oy=±l 



Vrs = _ [1 + { _ 1) M r+Pr ,- Prl] 



= 1 ^ gif (l-<Tr)(M r 



+Prt~Prl) 



(83) 



T r = ±l 



Using the projection operators, the partition function 
can be written as 



Z = Tr[e-^ H P C V S ] . 



(84) 



A Euclidean path integral representation can be obtained 
by splitting the exponential into M number of time slices. 



Z = Tr[{e-^V c V s ) M ] 



(85) 



where e = /3/M. Now the partition function can be writ- 
ten as 



/oo oo 
Y[df ia df ia d(j) ic d(j)i s E E E E e_i 

ia Ni = — oo Mi = — oo (Ti=±l Ti=±l 

(86) 

Here i = (r, r) runs over the 2+1 dimensional space-time 
lattice with r = 1,2, ...,M time slices. The action, 5, 
has the following form 

M 

S = S* + Sf" +Sp +eJ2 H ( N r, M T , <p 

r=l 

(87) 



with 



M 



r.T—1 a 
M 

S*° = 53 JV T [0 TC - T _i, c + |(1 - r T )] , 

r.r— 1 

M 

r.T—1 

(88)' 

Here the spatial index r is suppressed for clarity. The 
Ising variables a T and r T are defined on the links con- 
necting adjacent time slices and can be regarded as the 
time component of the Ising gauge fields. 

The sum of H t and Ha, can be decoupled using the 
Hubbard-Stratanovich fields Xrr' and rj rr i. 



-c(b*+b*) = J TJ Y[d X rr'dx; rl dr, rr ,dr,; r ,e 



* St,. 



Using the expressions for H t and Ha 

H t = -tJ2 (btcfir'af&Va + h.C.) , 
rr' ,a 

H A = A 53 arr'(bt a b r 'afrafr'a + h.C.) . (90) 



we have 



St ' A = \ e E { [ 2 IXrr'| 2 - Xrr'(Ka b r'a + tfrafr'c 



+ 



-h d rr > A.f ra f r 



2\f] rr '\ Tj rr > {b ra b r > a tfrafr' 



CL rr > A frafr' a) 

Rearranging terms, we get 



(91) 



>t,A 



2\ Xrr >\ 2 +2|r? rr /| 2 
-t(x rr' TJ rr ' ) fra. fi 



r a 
r'ct 



G>rr' ^(Xrr' ?]rr'^) frafr' a H~ CC. 

Rewriting this in terms of 6 r and z r , we get 



(92) 



2| X rr'| 2 + 2|r?„'| 2 



t(Xrr f Vrr'^frafr'a 

a rr ' A^Xrr' ^]rr') fra.fr' a 



~(Xrr> + T]rr')(blbr'z}.Z r ' + b\b r >Z r z\ t ) + C.C.j . (93) 

In order to decouple b r from z r , another Hubbard- 
Stratanovich transformation is necessary. Using similar 
procedures, the following term 



l e ^2i(Xrr' +11rr>)(b r br>Z r : Z r > + b\b r , Z r Z f r ,) + C.C. 



(94) 



2\\ rr ,\ 2 + 2\£„,\ 2 + 2\p r 



can be decoupled as 

- J^ e ^2(Xrr> +Vrr>) 
rr' 

-)~2|^ rr '| [\ rr r ~\- ^ rr ' ^Z r Z r ' {Prr' ~\~ Qrr'~)^r^ r ' 
i^Prr' Qvr' ^rr ^rr'^^^^r' C.C. 



(95) 



We now make a saddle point approximation and keep the 
Ising fluctuations around this saddle point. The natural 
choices are 



Xrr' T]rr' &rr'7~rr'Xf 5 
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(Xrr' + »7 rr ')(A rr .' + £rr' + p* r / + g* r ,) = OVr'Xs ■ (96) 

where oy r / = ±1 and T rr / = ±1 are Ising fluctuations. We 
drop all of the constant terms and define the following 
variables 

1 1 1 1 . . . 

*/ = 4 f X/, t c = — xc, t. = —Xs, t A = -A Xf (97) 



e v " 1 T ' = e^'-j 

Mi 

with 



S T B = -iN [2^ -£(1-7*. 



Here L T „- is defined as 



to obtain 



$ c 1 

*j , _ 

2?r 2 



(102) 



(103) 



(104) 



■err 



EE 



(Trr'Trr'itffraefr'a ~t~ ^A^rr'frafr'a) 



with $■ 



|r (1 — Tij) is the gauge invariant 



c'rr' 

bib r , + t s 



(98) 



Combining all the results, the approximate full parti- 
tion function can be written as 



oo oo 



E E n E E e_s ' (^) 

Ni=—oo M ( =— oc (ij) Uij=±l ry=±l 

where er^ and are Z 2 gauge fields living on the nearest 
neighbor links of the space-time lattice. The total action, 
S, is given by 



S = S f r + Sf° + Sf° + S A + So + S a + S v (100) 



with 
H 

si- 



% ^ , ^ ^ [fiajp'ij^'ij fja fia)] > 



■i J2 W 

i^j—i — r 

■i M < 



7T 



ic V]c , 2 



+ 77( 1 -^)] - 



7T 

'is - Pjs + ^(l - O'yjJ ! 



Sa — ( ^ ' tA&ij T ij{ a ij fiafju + C.C.) 



So 



\tf a ij T ij fia fja + ^c^jj &j &} 



i,j—i-\-x ot 



+ t 8 (JijZ*Zj + C.C. 



i 

S v = evY^(Mi) 2 , 



(101) 



where f and x represent the time and spatial links, dy = 
a rr ' on the spatial links and zero otherwise. 

Using the Poisson resummation formula, one can show 
that 



N, 



phase difference across the temporal link. Int denotes 
the integer part. One can see that the Berry phase term 
for (Tij is absent. This is due to the fact that we have 
equal amplitudes for up-up and down-down pairing in 
the equal spin pairing state, analogous to particle-hole 
symmetry in the charge sector. 

Gathering these terms, the final form of the action is 
given by 



S — Sf~\~Sc~\~Ss~\~Sj^~\~Sg 



(105) 



with 



Sf — ^ ] (J ij T ij\t{jfiafja + t Q-ijfiafja + C.C. 



^ ' fia fia ; 

ia 

S c = -Y / t c ij T ij (b*b j + c.c.) , 

ij 

-^tfjCTijiZiZj + C.C.) . 



S. 



(106) 



Here is et c on the spatial link and on the temporal 



link. Similarly t\j is et s on the spatial link and on the 



temporal link, t jj = etf on the spatial link and tfj = — 1 
on th e temporal link. And t A = etA- The last term of 
(105) corresponds to the Maxwell terms for the Z 2 gauge 
fields, that we assume are generated after we integrate 
out excitations at high energies. 



S g = - KiYJ] - K 2 YJ] 



□ □ 

□ □ 



□ □ 



(107) 



These are the simplest terms providing dynamics of the 
gauge fields that are consistent with the gauge symme- 
tries 



Z 2 t • b h ► t^bi, f ia ^ tifia? ^ ^itjT^ 

Z 2 a • %i ^ SiZii fia * ^ifia: &ij * SiSj&i- 



(108) 



where tj and Sj are ±1. In the future we will call any 
particle that transforms under the first and the second 



transformations of (108) as having Z 2r and Z 2a charges 
respectively. 
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B. Spin Singlet Insulating Phases 

Before discussing possible spin singlet insulating 
phases of the combined action ( |l05|) - (107) it is use- 
ful to review properties of a pure Z2 X Zi gauge theory 
(107). Under duality transformation defined in |7|, |44| , |45[ 
this model becomes a generalized Ashkin- Teller model 

Pel 



Sat = - Kdi ^ ViV i ~ Kd2 ' 
(v) (*i> 

- K d3 2J UiViUjVj 

m 



(109) 



Here itj and Vi are Ising variables defined on the dual 
lattice in d = 2 + 1. We can identify five possible phases 



of (|109| ): 



1) Fully ordered phase (u) 7^ 0, (v) 7^ 0, (™) ^ 0; 

2) Partially ordered phase (u) 7^ 0, (w) = 0, (uv) = 

3) Partially ordered phase (u) — 0, (v) 7^ 0, (uv) = 

4) Partially ordered phase (u) — 0, (v) — 0, (m>) 7^ 

5) Disordered phase (u) = 0, (v) = 0, (uv) = 0. 

As pointed out in [0 the Ising variables of ( |109| ) cor- 
respond to the vortices of the original gauge model 
They describe gauge field configurations with plaquette 
products equal to minus one, i.e. plaquettes pierced by 
Z2 fluxes. Following |t]] we call such Z2 vortices "visons". 
In fact we have three kinds of visons: cr-visons that de- 
scribe Z2 vortices of a, t- visons that describe Z2 vortices 
of r, and [crr]-visons that describe a composite of a and 
t Z2 vortices. The three are not independent, any one 
of them can be thought of as a composite object of the 
other two. However, we should treat all of them on equal 
footing since they represent distinct topological objects. 
The appearance of the long range order in the Ashkin- 
Teller model corresponds to the condensation of visons in 
the original gauge model and describes transition to the 
confining phase. From these arguments it follows that 
t here are five distinct phases of the pure gauge model in 
(107). One fully confining phase, three partially confining 
phases, and one fully deconfining phase, that correspond 
to the fully ordered, three partially ordered, and one fully 
disordered phases of the Ashkin- Teller model. 

1) Fully confining phase, a- and r- visons are condensed 
simul 

taneously. This also implies condensation of [ut]- 
visons. 

2) Partially confining phase, r-visons are condensed and 
<7- and [<tt] visons are gapped. 

3) Partially confining phase, a visons are condensed and 
r and [<rr]-visons are gapped. 



4) Partially confining phase, [err] visons are condensed 
and a and r visons are gapped. 

5) Deconfining phase. All visons are gapped 

Condensation of visons has dramatic effects on the mo- 
tion of spinons, holons, and neutral fermions in the model 
( |105| )-( |l07"| ). We find drastically different excitation spec- 
tra depending on what vortices are condensed. The rea- 
son for this is a geometrical phase factor of it that parti- 
cles with i^-charges acquire when they circle around an 
appropriate Z2 vortex. For example, spinons and neutral 
fermions get a geometrical phase factor of n when they 
are transported around a cr-vison, and holons and neu- 
tral fermions get a minus sign when they circle around 
a r-vison. This means that when visons are present in 
the ground state, the coherent motion of the correspond- 
ing particles is highly frustrated and they may not be 
considered as elementary excitations. Only the particles 
that are neutral with respect to the appropriate Z2 sym- 
metry may propagate freely in a phase with condensed 
visons. And the particles that carry such Z2 charges will 
have to bind into neutral pairs. This is the essence of the 
confinement argument discussed in fl4|,|7||. 

When we apply the geometrical phase - confinement 
argument to the spin singlet insulating states we find the 
same phases as discussed in Section [V E| . 



• In a phase of type 1) all kinds of Z2- vortices are 
condensed. Therefore, particles that carry any Z2 
charges will be bound. This is a fully confining 
phase where only fully neutral composites are al- 
lowed. Holons, spinons and neutral fermions are 
confined. Phase CSf. 

• In a phase of type 2) we have a condensate of r- 
visons. As a result particles that carry Z2 T charges 
are confined, but particles that carry Z20 charges 
are liberated. Spinons are free, and holons are 
bound to the neutral fermions. Phase CfSb. 

• In a phase of type 3), that has a condensate of cr- 
visons, we have a confinement of particles with Zi a 
charges and deconfinement of particles with Z2 T 
charges. Holons are free, and spinons are bound 
to the neutral fermions. Phase CbSf. 

• In a phase of type 4) we do not have individual a- 
and r- visons in the ground states, but only their 
composites, [err] -visons. The geometrical phase ar- 
gument becomes somewhat subtle when we con- 
sider [cr] -visons. Particles that carry either one of 
Z2T or Z2 a charges will get a ir phase shift when 
they circle around such a vortex. However, par- 
ticles that carry both charges acquire no phase. 
So, in a D-type phase particles that carry one of 
the Z2T or Z2a charges are confined, but particles 
that carry both charges are deconfined. Holons and 
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spinons are bound, and neutral fermions are free. 
Phase CSbNf. 

• Finally in a phase of type 5) we have no condensed 
visons, which means that all the particles are lib- 
erated. Holons, spinons, and neutral fermions are 
deconfmed. Phase CbSbNf. 



C. Broken Symmetry Phases 

In this section we show using Z2 x Z2 theor y th at eve n 
states with the long range order in the model (105)- (107), 
i.e. p wave superconductors, spin singlet superconduc- 
tors, nematic insulators, and nematic superconductors 
may differ in their topological ordering and carry the 
remnants of the spin-charge separation that appears so 
dramatically in the insulating phase. 

We begin by reviewing the case of a p wave supercon- 
ductor. 



• The simplest p wave superconductor that may be 
deduced from the model (105) - (107) is when 



holons and spinons condense simultaneously, so the 
system acquires hnite expectation values of b and 
z. The geometrical phase argument when applied 
to this system tells us that an isolated hc/2e vortex 
or a meron are no longer well defined excitations, 
since they acquire a phase shift of ir when circling 
around a holon or a spinon respectively. However, 
if we bind an hc/2e vortex with a r-vison we find 
that this composite can propagate freely. The geo- 
metrical phases acquired by the two upon encircling 
a holon add up to or 2ir. Equivalently, a meron, 
when bound to a cr-vison, becomes a well defined 
excitation in the presence of spinon condensate. 

• Another possible phase of a p-wave superconductor 
is when we condense holon pairs and spinons, i.e. 
b 2 and z. In this phase merons are still bound to 
cr-visons, however hc/2e vortices and r-visons are 
now deconfmed. The original holons are reduced to 
Ising variables, which we can call 6-isons, following 
[fjj . They carry the leftover of the charge symme- 
try, that was broken from U(X) to Z2, and are well 
defined excitations in this phase. 

• Analogously to the previous case we can consider a 
situation with condensed b and z 2 . This phase will 
have bound hc/2e vortices and r-visons and liber- 
ated merons and cr-visons. Spinons become Ising 
variables, z-isons, that carry the residual Z2 spin 
quantum numbers 

• A different type of a p-wave superconductor occurs 
when holon pairs and holon-spinon composites con- 
dense simultaneously, i.e. b 2 and bz acquire expec- 
tation values (this also fixes the expectation value 



for z 2 ). In such a phase spinons and holons are 
reduced to a single Ising variable, since knowing 
b automatically gives z. This 6z-ison carries the 
residual spin-charge quantum number of the sys- 
tem. A stable topological object in this phase may 
be constructed by taking any two of the set (hc/2e 
vortex, meron, r-vison, cr-vison). 

• Finally, we may have a phase with condensed holon 
and spinon pairs, b 2 and z 2 . This gives us separate 
6-isons, z-isons, hc/2e vortices, merons, r-visons, 
and cr-visons. 

The last four phases are the triplet analogues of the exotic 
SC* phase discussed in fjj in the case of singlet super- 
conductors. 

We now consider the case of a spin-singlet supercon- 
ductor. 

• The simplest kind of a spin singlet superconduc- 
tor occurs when we condense simultaneously holons 
b and cr-visons. The former ensures confinement 
of hc/2e vortices and r-visons, whereas the latter 
gives rise to binding of neutral fermions to spinons. 

• Another possibility is to have a condensate of holon 
pairs b 2 and cr-visons. This liberates hc/2e vortices 
and r-visons, produces 6-isons that carry charge 
Z2 number, and leaves neutral fermions bound to 
spinons. 

• Another option is to have a condensate of bosons 
b with gapped cr-visons. This means bound 
hc/2e vortices and r-visons, and liberated neutral 
fermions and spinons. 

• The most intriguing phase in this series is obtained 
when we condense holon pairs b 2 and holon-cr-vison 
composites. Excitations in this phase will be any 
pair from the set (hc/2e vortex, r-vison, ip, z) and 
6-isons. 

• Finally we can have a condensate of b 2 and gapped 
a visons. This gives unconfined hc/2e vortices, r- 
visons, neutral fermions, spinons, and 6-isons. 

Of the five phases above, four of the last ones may be 
considered as SC* phas 
es. 

The construction given above for p-wave superconduct- 
ing states and spin singlet superconducting states may 
be generalized to the case of spin-nematic insulators and 
nematic superconductors. In those cases, just as in the 
two discussed above, we find five possible states. One of 
these is a traditional version, whereas the other four are 
of the unconventional *-variety that may be thought of 
as containing traces of quantum number separation. 
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The reader may be worried that we do not find 
hc/Ae -vortices- 7r-disclinations in our discussion of var- 
ious phases of p-wave superconductors. As in the previ- 
ous sections we assumed that these excitations have been 
gapped out (see discussion after equation (|79|)). 

VI. DISTINGUISHING DIFFERENT 
FRACTIONALIZED PHASES 

In previous sections, we have seen how various frac- 
tionalized phases can arise in the context of Kondo lat- 
tice models and systems with a tendency towards p-wave 
superconductivity or superconductivity coexisting with 
magnetism. These phases can be described in the lan- 
guage of vortex and skyrmion condensation or in terms 
of a Z2 x Z2 gauge theory. However, one might wonder 
if these results are an artifact of these formalisms. In 
particular, one can ask how these phases can be distin- 
guished - both as a matter of principle and as a prac- 
tical experimental issue - from each other and from un- 
fractionalized phases. As Wen jll| and, more recently, 
Senthil and Fisher Q have emphasized recently, their 
'topological order' - i.e. the sensitivity of the ground 
state to changes of the topology of the system - provides 
one means of distinguishing fractionalized phases. 

This characterization of fractionalized phases is cru- 
cial because other heuristic definitions of fractional- 
ized phases can fail. To see why this is so, consider 
the intuitively-appealing statement that a fractionalized 
phase is distinguished from a conventional phase by ask- 
ing for the lowest energy excitation with, for instance, 
spin-1/2. In the conventional case, this would be an elec- 
tron which also carries an electric charge e. In the frac- 
tionalized phases of the kind discussed above, one might 
expect that the corresponding excitation is a spinon 
which is charge neutral. However, this test for fractional- 
ization fails if there is an attractive interaction between 
the holons and spinons which binds them into an electron 
at low energies. This could, in principle, happen without 
going through a phase transition. (Unlike in an unfrac- 
tionalized phase, holons and spinons would still exist as 
unbound excitations, but at higher energies.) Then, the 
lowest energy excitation with spin-1/2 is an electron (as 
opposed to a spinon) though the system is adiabatically 
connected to a fractionalized phase (see Ref. [^3| , p8|j29| l for 
a discussion of this effect). Furthermore, other tests such 
as the vanishing of the quasiparticle residue at some point 
of the Brillouin zone also fail in this situation. Hence, we 
turn to the characterization in terms of the topological 
properties of the system. 

Topologically ordered systems are partially character- 
ized by their ground state degeneracy on the annulus, 
the torus, or higher-genus manifolds, over and above 
any degeneracy which may be due to broken symmetry. 
Consider the CbSf phase. It has a two-fold degenerate 



ground state on the annulus. The two ground states cor- 
respond to periodic and antiperiodic boundary conditions 
for holons and spinons as they encircle the center of the 
annulus. In either case, electrons themselves have peri- 
odic boundary conditions, as they must. In an unfrac- 
tionalized phase, spinons and holons are confined within 
an electron so the two states are identical; the excita- 
tions which could distinguish them are not part of the 
spectrum. By the same reasoning, the CfSb and CSbNf 
phases also have two degenerate ground states on the 
annulus. By extension, all of these states have ground 
state degeneracy 4 9 on a genus g surface. On the other 
hand, CbSbNf has four degenerate ground states. We can 
independently choose periodic or antiperiodic boundary 
conditions for the charge and spin bosons. The bound- 
ary conditions for the neutral fermions are then deter- 
mined by the requirement that electrons must have pe- 
riodic boundary conditions. On a genus g surface, it has 
degeneracy 16 9 . 

These degeneracies can be interpreted in terms of the 
vison spectra of the fractionalized states. The two ground 
states of CbSf on an annulus correspond to the presence 
or absence of a r vison (i.e. a v) in the center of the an- 
nulus; the two ground states of CfSb correspond to the 
presence or absence of a a vison (a v'); the two ground 
states of CSbNf correspond to the presence or absence of 
a err vison in the center of the annulus. The four ground 
states of CbSbNf correspond to the presence or absence of 
a and r visons in the center of the annulus. The interpre- 
tation of these ground states in terms of visons forms the 
basis for an experimental probe of topological order pro- 
posed by Senthil and Fisher |l9fl . We will return to this 
issue later but let us, in the meantime, continue to pur- 
sue the question of the distinction in principle between 
different fractionalized phases. 

Different states at the same level of fractionalization 
have the same ground state degeneracy; CbSf, CfSb, and 
CSbNf all have two degenerate ground states on the an- 
nulus. In order to distinguish them, we must consider 
their quantum number spectra. CSbNf does not have 
spin-charge separation, i.e. it is not possible to isolate a 
charge-0, spin-1/2 excitation at finite energy cost. Fur- 
thermore, it is possible to isolate a neutral fermionic exci- 
tation. Both of these stand in contrast to CbSf and CfSb 
which exhibit spin-charge separation but do not support 
neutral fermionic excitations. Hence, we conclude that 
CSbNf is distinct from the other two states despite hav- 
ing the same ground state degeneracy. 

One might be tempted to conclude that CbSf and CfSb 
are distinct because the lowest-energy charged excitation 
is a boson in one phase and a fermion in another phase. 
However, if a holon in CbSf forms a bound state with 
a t vison, the resulting bound state will be fermionic; 
similarly, if a spinon in CbSf forms a bound state with a 
r vison, the resulting bound state will be bosonic. Hence, 
as a result of the seemingly innocuous formation of bound 
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states, the CbSf and CfSb states appear to morph into 
each other. Thus one is instead tempted to conclude that 
the CbSf and CfSb phases can be adiabatically connected 
to each other. 

This contention is supported by considering the singlet 
superconducting state which results if holons condense 
in CbSf or if holon-er vison composites condense in CfSb 
(see Figure . It is easy to see that the superconducting 
states in either case are conventional and are smoothly 
connected to a BCS state. The superconducting state can 
be disordered by vortex condensation. This will yield a 
fractionalized state (with a two-fold degenerate ground 
state on an annulus) if vortex pairs condense but indi- 
vidual vortices are uncondensed. Since the result could 
be either CbSf or CfSb, this appears to support the pos- 
sibility that there is no phase boundary between these 
phases in the part of the phase diagram near the singlet 
superconducting phase. 

However, there is a logically-possible alternative, 
namely that an operator which is irrelevant in the super- 
conducting phase and at the critical point becomes rele- 
vant at the fixed points characterizing the fractionalized 
phases. In that case, the actual nature of the resulting 
fractionalized phase depends on short distance physics 
- the value of the coupling which is formally irrelevant 
in the superconductor - and is not uniquely dictated by 
knowing that there is proliferation of hc/e and with hc/2e 
vortices remaining gapped. 

Despite this caveat, a scenario in which CbSf and CfSb 
are smoothly connected to each other in the vicinity of 
their transition to the superconducting state is appeal- 
ing and plausible. This does not necessarily mean that 
CbSf and CfSb are not distinct phases. Their relationship 
could be similar to that between a liquid and a gas, which 
are separated by a first-order phase transition line which 
terminates at a critical point, beyond which a liquid and 
a gas can be adiabatically connected without crossing 
a phase transition line. In appendix C, we show that 
precisely such a scenario does occur in simpler (though 
somewhat different) Z2 x Z2 gauge theory models. Thus, 
we tentatively suggest that the first-order phase transi- 
tion between the CbSf and CfSb phases terminates at 
a critical point. Beyond this critical point, there is no 
distinction between these phases, and it is in this region 
of the phase diagram that there is a phase transition to 
the superconducting phase. 

VII. FLUX-TRAPPING EXPERIMENTS 

Let us now consider the practical issue of how we can 
identify whether a given system in an unknown phase is 
fractionalized or not and, if it is fractionalized, then what 
its fractionalization pattern is. 

To proceed, note first that the CbSf phase contains 
in it the seed of superconductivity. As argued in Ref. 



jjj, condensing the charged boson provides a natural 
non-pairing route to superconductivity (of a conventional 
kind). Similarly, the CfSb phase contains in it the seed 
of magnetism - simply condensing the spinon leads to a 
conventional state with some kind of magnetic long range 
order. However, it is possible to imagine a transition be- 
tween the CfSb phase and a superconductor which occurs 
when a composite formed by a holon and a a vison con- 
denses. Similarly, it is possible to imagine a transition 
between the CbSf phase and a magnetic phase which oc- 
curs when a composite formed by a holon and a a vison 
condenses. 

The feature of most interest to the following discussion 
is simply that a direct phase transition should be possible 
between the CbSf and CfSb phases and a conventional 
superconductor. Upon going through such a phase tran- 
sition, the visons of these phases acquire hc/2e units of 
electromagnetic flux to become the hc/2e vortices of the 
superconductor. This may be exploited to devise a sen- 
sitive test for the topological order in the CbSf phase, 
as argued in Ref. Jl^Q . 

The test proceeds as follows. Consider an annular sam- 
ple of a material which is in a conventional superconduct- 
ing phase and let us suppose that we can tune the sample 
parameters adiabatically so that the sample makes tran- 
sitions between the superconducting phase and the CbSf 
and CfSb phases. Suppose that hc/2e units of electro- 
magnetic flux are trapped in the annulus when the system 
is in its superconducting phase. There must also be a vi- 
son trapped in the annulus so that the holon condensate 
can have periodic boundary conditions (without which it 
would cost infinite energy) : the antiperiodicity caused by 
the flux hc/2e is cancelled by the antiperiodicity due to 
the vison. If the system is taken into the CbSf phase, the 
flux escapes since there is no holon condensate trapping 
it, but a vison will remain since it will cost energy (the 
vison gap) to unwind the antiperiodic boundary condi- 
tions of the (neutral) spinons. If the system is returned 
to the superconducting state, then it must generate flux 
±hc/2e so that the holon condensate can again have peri- 
odic boundary conditions. The same analysis holds if we 
take the system into the CfSb phase except that we have 
to replace 'holon' in the above description by 'holon- vison 
composite'. On the other hand, if the system undergoes 
a transition to an unfractionalized phase, then the vi- 
son can escape since there are no deconfined spinons or 
holons whose boundary conditions would be affected by 
its escape. 

Of course, this experiment would simply be confirm- 
ing the result which we arrived at in the previous section: 
that the CbSf and CfSb phases can be adiabatically con- 
tinued into each other, particularly in the neighborhood 
of a conventional singlet superconducting phase. 

Let's now consider a more complicated flux-trapping 
experiment in which, as an intermediate step, we take 
the system through the higher-level fractionalized phase, 
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CbSbNf (see Figure [j]). This phase has two distinct vison 
excitations. One of these visons can be envisioned as a 
descendent of the r vison of the CbSf phase; we will refer 
to this as v. The other can be envisioned as a descendent 
of the a vison of the CfSb phase or as a by-product 
of the further fractionalization of the fermionic spinon 
of CbSf; we will refer to this as v' . A direct transition 
from CbSbNf to the CbSf phase occurs when the visons 
v' condense while that from CbSbNf to CfSb occurs 
when the visons v condense. The presence of two distinct 
visons in the CbSbNf phase distinguishes it from the 
CbSf and CfSb phases - indeed it will have a ground 
state degeneracy of 16 on a torus. 

Now consider a conventional BCS superconductor. 
This is obtained from CbSf by condensing the holon. 
The flux-trapping experiment performed by moving be- 
tween the superconductor and CbSf gives a positive re- 
sult. Now consider a modification of the experiment so 
that we start in the superconducting phase, move first 
to CbSf, then to CbSbNf, then back into CbSf before 
finally going back into the superconductor. This again 
gives a positive result since v is trapped in the annu- 
lus and it can never escape. Upon making the transi- 
tion between the CbSf and CbSbNf phases, a v' will be 
generated with probability 1/2 since the ground state of 
CbSf with one v will make a transition to either of the 
corresponding ground states of CbSbNf with equal prob- 
ability. However, this v' will escape upon the transition 
from CbSbNf back to CbSf. Now consider a further 
modification in which we go all the way from the su- 
perconductor to the CfSb phase through the CbSf and 
CbSbNf phases and then return by the same route to 
the superconductor. The result of this experiment will 
be negative half of the time. This is because in going 
from CbSbNf to CfSb, the vison v condenses. Thus v 
which was trapped in the hole until the phase CbSbNf 
was reached can escape on moving into the CfSb phase. 
In going from CfSb back to CbSbNf, a v is generated 
with probability 1/2 - the two ground states are obtained 
with equal proability. This v, if it is generated, will lead 
to the generation of flux hc/2e in the superconducting 
state. 

Hence, there appears to be a difference between the 
CbSf and CfSb which can be detected in this exper- 
iment. It appears that these phases cannot be contin- 
uously connected - since the probability of a negative 
result for the flux-trapping experiment of the previous 
paragraph must jump from to 1/2 - at least in the 
vicinity of the CbSbNf phase. This can be understood 
in the following terms. In the CbSbNf phase, there are 
two distinct types of visons, v and v' . If one or the other 
condensed, a transition occurs to CfSb or CbSf. The re- 
maining vison in CbSf 'remembers' that it is a v vison. 
Meanwhile the vison in CfSb 'remembers' that it is a v' 
vison. However, if we take the system far from CbSbNf 
so that a bound state can form between a v and a holon 



and also between a v and a spinon, then v now looks 
like a v 1 and the distinction between the two phases is 
blurred. Combining this reasoning with that of the pre- 
vious section, we propose the phase diagram of figure @. 




FIG. 6. A schematic phase diagram indicating how the 
CbSf, CfSb, CSbNf, CbSbNf, and conventional singlet su- 
perconducting phases might fit together. The thick lines are 
first-order phase transition and the thin lines are second-order 
phase transitions. 

VIII. DISCUSSION 

When electrons interact strongly, a number of inter- 
esting phenomena are known to occur, including uncon- 
ventional superconductivity and magnetism. As we have 
seen in this paper, many of the physical settings which 
give rise to these phenomena also have the potential to 
exhibit electron fractionalization. Different theoretical 
approaches, adapted to these specific systems, suggest 
seemingly different fractionalized phases. It is natural to 
ask if these phases are truly different and, if so, what 
their organizing principle is. 

In this paper, we have pursued the idea that a 
crisp and coherent way of understanding quantum num- 
ber fractionalization is provided by the concept of topo- 
logical order introduced in the context of fractional quan- 
tum Hall liquids pC|| , anyon superfluids |4S}| , chiral spin 
states p9[ , and short range resonating valence bond spin 
states ]50|] . We presented two approaches for understand- 
ing such topological order. The first one relies on the 
recently developed Z2 gauge theory of spin-charge sepa- 
ration, originally suggested for the high T c cuprates, and 
generalizes it to a Z2 x Z2 theory to include possible 
fractionalization of spin and charge quantum numbers. 
Some of the interesting fractionalized phases are: CbSf 
(bosonic holons and fermionic spinons), CfSb (fermionic 
holons and bosonic spinons), CSbNf (bound bosonic 
holons and spinons and neutral fermions), and CbSbNf 
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(bosonic holons and spinons and neutral fermions). Any 
one of these phases can be further characterized by pos- 
sible broken symmetries with conventional order param- 
eters. Each of the fractionalized phases corresponds to 
a different dcconfining phase of the pure Z 2 x Z 2 gauge 
theory and will have appropriate topological Z 2 vortices, 
visons, as finite energy excitations. 

An alternative picture of fractionalization which is also 
presented in this paper uses the language of quantum 
disordered superconductors and magnets. When topo- 
logical ordering - defined by the suppression of certain 
defects - occurs, the Goldstone modes associated with 
various broken symmetries can screen the corresponding 
quantum numbers of the fermionic quasiparticles. In this 
way, these quasiparticles can be bleached of some or all 
of their quantum numbers. This may be implemented 
mathematically with U{\) particle- vortex duality in both 
the charge and spin sectors. We arrive at essentially the 
same picture as that of the Z 2 x Z 2 gauge theory. In 
those insulating phases in which hc/2e vortices are con- 
densed, charge is bound to the fermionic quasiparticles. 
When hc/2e vortices are gapped and hc/e vortices are 
condensed, charge carrying "holons" can propagate sep- 
arately from the electrically neutral fermionic quasipar- 
ticles. In the spin sector, we can consider either meron 
or skyrmion (which carry twice the topological charge of 
merons) condensation, with gapped merons in the latter 
case. In the former case, spin is confined to the fermionic 
quasiparticles, and in the latter case "spinons" will ex- 
ist as independent objects, deconfincd from the fermionic 
quasiparticles. We have also discussed the possibility of 
quantum disordered phases in which the condensed topo- 
logical objects are the hc/2e vortex - meron composites, 
but not hc/2e vortices or merons separately. Such phases 
have "spinons" and "holons" bound together but decon- 
fined from the neutral fermionic quasiparticles. 

An important issue discussed in this paper is whether 
one can distinguish the phases obtained by quantum 
disordering the spin and charge sectors of the system, 
for example the phases CfSb and CbSf of the quantum 
disordered p-wave superconductor. The simplest choice 
seems to be the identification of the spin excitation as 
a fermionic or bosonic particle. This, however, is not a 
reliable tool. In the Z 2 x Z 2 gauge theory formulation, 
both spinons and holons carry Z 2 charges, so a bound 
state of a Zi vortices with either one of them (this can 
also be thought of as attaching Wilson loops to the par- 
ticles) will change its statistics from fermionic to bosonic 
or vice versa [ p^||39|j40| , |5l| | . In the deconfining phase such 
vortices are gapped. However, if a bound state between 
a Z 2 -charge carrying particle and a Z 2 vortex forms, this 
bound state may have a lower energy than the original 
particle. This means that in both CfSb and CbSf phases 
the lowest energy spin- or charge-carrying excitations can 
exist as either bosons or fermions. The subtleties dis- 
cussed above lead us to consider flux-trapping experi- 



ments of the type discussed in Section VII. Combining 
all of these considerations, we outlined one scenario in 
which CbSf and CfSb phases can be separated by a first 
order transition which terminates at a critical point. On 
the other hand, one can go from CbSf to CfSb through 
CbSbNf phase by two continuous transitions. Thus, if this 
scenario is correct, the relation between CbSf and CfSb is 
somewhat similar to that between liquid and gas phases. 
We however defer offering any definitive conclusion. 

Spin charge separation in one dimensional systems is 
fundamentally different from its two dimensional coun- 
terpart, since it does not involve topological order. An- 
other nontrivial realization of electron number fraction- 
alization which is analogous to that presented here can 
occur in multi-component quantum Hall systems and was 
discussed in references 37 5j|. 



Another avenue for further research is the investiga- 
tion of quantum-disordered states of triplet supercon- 
ductors with more complicated spin structures appear- 
ing in some of the superfluid phases of 3 He. We expect 
that these will share some features of non-collinear spin- 
density- waves |52). Further exotic phases are likely to 
occur upon quantum-disordering states with multiple or- 
der parameters. We have considered one of the simplest 
cases of this - antiferromagnetism and superconductivity 
- but there are more complicated possibilities, involving 
incommensurate charge and/or spin order. 

In addition to the phases CbSf and CfSb that have ap- 
peared in the literature bef ore, we proposed the possibil- 
ity of two additional quantum number separated phases 
in these systems: phase CSbNf in which the excitations 
are a spin-1/2, charge e boson, a neutral spinless fer mion 
and a vison and phase CbSbNf with a charge e spin- 
less boson, a neutral spin 1/2 boson, a neutral spinless 
fermion, and two distinct visons. simultaneously con- 
densed hc/e vortices and skyrmions. 

The possibility of a higher SO (5) symmetry which 
unifies d-wave superconductivity and antiferromagnetism 
has been suggested for the high T c cuprates and organic 
superconductors in f54j| . In the SV 2 i?ii04 materials, a 
similar symmetry has been proposed in [pa] which com- 
bines p-wave superconductivity and ferromagnetism. An 
effective model for the coupling of quasiparticles to a 
fluctuating SO (5) order parameter has been derived in 
|56f . In this model holons and spinons are not segregated 
into independent quasiparticles from the very beginning 
but are naturally combined into composite quasiparti- 
cles which transform as spinors of SO (5). Such spinors 
are spin doublets and carry charge e |57],|58|]. There are 
also neutral fermions which carry no quantum numbers. 
One can see a striking resemblance between these exci- 
tations and the excitations in the phase CbSbNf. This 
suggests the interesting possibility that the restoration 
of the SO(5) symmetry in models with strong quantum 
fluctuations manifests itself not in the existence of a bi- 
critical point on the phase diagram, but in the appear- 
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ance of a specific form of quantum number separation of 
the electrons. A detailed discussion of quantum disorder- 
ing phenomena in models with SO (5) symmetry requires 
a detailed analysis of the non-Abelian Berry's phases in- 
volved in the description of SO(5) spinors and will be 
presented in subsequent publications. 

The states which we have discussed, as well as the more 
complicated ones alluded to above, have potential appli- 
cation to a variety of materials, including not only the 
cuprates p2| |, but also Sr 2 RuO i heavy fermion su- 
perconductors, such as Celn^ p4|; and organic super- 
conductors, such as k — (ET) 2 Cu[N(CN) 2 ]Cl |§. All 
of these compounds have magnetic (in come cases incom- 
mensurate) phases in proximity to p-wave or d-wave su- 
perconducting states. It is possible that pressure, chem- 
ical substitution, magnetic field, etc. might drive a tran- 
sition into one of the phases described here in which the 
magnetism and the superconductivity are disordered by 
quantum fluctuations. 

Ideas presented in this paper should also apply to 
Bose-Einstein condensates of spinor bosons, such as al- 
kali atoms 23 Na and 87 Ru which have a hyperfine spin 
F = 1. For example, when restricted dimensionality or 
quantum fluctuations destroy the spin ordering we ex- 
pect to find condensation of pairs of atoms into a global 
spin singlet state, and when quantum fluctuations in the 
charge sector destroy the U(l) phase ordering we can 
find states characterized by a spin nematic order. Some 
of these phenomena have been discussed in | p5| . 

To summarize, we have studied the possibility of frac- 
tionalization in systems with ordering tendencies in the 
charge and spin sectors, including Kondo lattices, p-wave 
superconductors and systems with simultaneous d-wave 
superconducting and antiferromagnetic fluctuations. In 
the case of p-wave superconductors we find that the rich 
internal structure of their order parameter allows for the 
existence of the following quantum disordered phases: 
a charge 4e singlet superconductor, a spin singlet in- 
sulator, and a spin nematic insulator. For both the 
p wave superconductors and the cJ-wave superconduc- 
tor/antiferromagnet systems, we find that the quantum 
disordered phases may have separated quantum numbers, 
depending on the topological order, which can be char- 
acterized by specifying the nature of the finite energy Z 2 
visons. 
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APPENDIX A: KONDO LATTICE MODEL 

In this appendix, we provide some of the details of the 
Z 2 gauge theory reformulation of the Kondo lattice model 
discussed in Section [□]. Consider the Hamiltonian in Eqn. 
[j]. As in the discussion of the pure exchange Hamiltonian, 
we first replace the spin operator S~ by the boson opera- 
tor b sr ~ e llpT . The exchange Hamiltonian takes the form 
Eqn. U and the Kondo coupling takes the form Eqn. |24|. 
The electron hopping term is unaffected. We now change 
variables to spinon and holon operators as in Eqns. [?], 
and p7| . The terms H t , and H ex are now given by 
Eqns. |29|, ^a, and ^ respectively. In the presence of the 
Kondo coupling between the local moments and the con- 
duction electrons, the total (z component of the) spin at 
each site is 

n r + -cl<7 z c r (Al) 

We therefore define the total spinon number 

N* ot = 2n r + c\a z c r (A2) 

Note that N* ot is conjugate to the phase cf> r of the spinon 
field. We will work with the operators (z r , N* ot , n^ r , ryx. r ) 
instead of the original electron and local spin S r opera- 
tors. This change of variables however introduces some 
redundancy - the Hilbert space of states on which the 
holon and spinon fields operate is larger than the physi- 
cal set of states. This may be seen by noting that with 
the definition above, the operator A^* ' must satisfy 

AT*°' - 4<j z c r = even (A3) 

^From the definition of the holons, it follows that 
c\(j z c r = n\a z r\ r . Furthermore nj,a z rj r has the same par- 
ity as r\\r] r . Thus we have the constraint 

N* ot - T]lr] r = even (A4) 

The Hamiltonian needs to be supplemented with this 
constraint to correctly represent the original model (be- 
fore the change of variables). 
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It is useful to rewrite the exchange and Kondo parts 
of the Hamiltonian as follows: 



H K + H ex = Jk^2 \ r lt] r lri + h. 



U 



<rr'> 



r r 

The last term is an interaction between the holons. 
Clearly this term cannot affect issues of confinement of 
the holons with the spinons. We will therefore drop it 
for the present discussion. The last but one term repre- 
sents an interaction between the spinon density and the 
holons. We again expect that such an interaction is also 
unimportant for issues of the stability of fractionalized 
phases. We will therefore drop this too. 

We may now derive a functional integral representation 
of the system, proceeding as in Ref. §. The resulting 
action is 



S — S f ~\~ S j- 



(A6) 



Here S T represents terms involving coupling along the 
(imaginary) time direction. This and the Berry phase 
Sb are exactly the same as in Ref. Rj. The spatial part 
of the action is 

S r = Si + S K + S n (A7) 
Si = -e ^ trr> [z\z rl + (vl^r) r i + r^jVlJ + h.c)j 

<rr'> 



(A8) 



S K = +cJk ^2 yiU T iri + c - c 

T 

Su = —ej 22 ( z r^ z r + ^.C. 



<rr'> 



We now combine the terms «S7 and Su and rewrite them 
as 



zj ^2 ( Z t Z r' + 



<rr'> 

+h.c. + 0( V 4 ) 



2J 



(A9) 



The last term is a four-holon interaction which we will 
ignore on the grounds that it cannot affect issues of frac- 
tionalization. It is convenient to further rewrite the ex- 
pression above as follows: 



Si +Su = — ^2 

<rr'> 



-,' + -Jj'l'r'l' 



h.c. 



u rr ' J. 

2J*' 



h.c 



(A10) 



We may now decouple each of these two terms with a real 
Hubbard-Stratanovich field to write 



-(Si+Su) 



[V X Vp]e 



<s x +s„) 



(All) 



S X = T E Xrr' ~ 2Xrr ' \ Z r Z r' + ^V^rVr' + k.C. ) (A12) 



<rr'> 



Sf> = ~7T ^ Plr' - 2 Prr> ( z\ Z r> ~ TTf^r' ~ h.C. ) (A13) 



<rr'> 



2J 



Note that Xrr' = Xr'r while p rr r — —p r i r . We now con- 
sider evaluating the x> P integrals in a saddle point ap- 
proximation. Looking for uniform saddle points, we write 



< Xrr' >= Xo; < Prr' >= PQ 



(A14) 



Note that a non-zero value of po requires specifying di- 
rections for all the links of the lattice. The saddle point 
equations are 

XO = < z t z r' + -^jrfrVr' + h - C - > ( A15 ) 

Po = < z\z T > - -^j-vhlr' - h.c. > (A16) 

Note that po must be pure imaginary as it is the expecta- 
tion value of an antihcrmitian operator. With non-zero 
po, the saddle point action therefore becomes complex - 
this breaks time-reversal symmetry (and possibly various 
lattice symmetries due to the need to specify directions 
to the links). We restrict ourselves to time- reversal in- 
variant saddle point solutions, and therefore set po = 0. 
The resulting saddle point action then preserves all the 
global symmetries of the original model. However, it does 
break the local Z2 symmetry introduced by the change 
of variables to the holons and the spinons. This can be 
remedied by keeping a particular set of fluctuations about 
the saddle point, namely those associated with a change 
in the sign of the fields Xrr' ■ 



Xr 



X0°r 



(A17) 



with o~ rr i = ±1. The er rr / may be identified as the spatial 
components of a Z2 gauge field. We thus finally arrive at 
the action in Eqn. 30 



APPENDIX B: QUANTUM NUMBER 
SEPARATION IN SYSTEMS WITH D-WAVE 

SUPERCONDUCTING AND 
ANTIFERROMAGNETIC FLUCTUATIONS 

Nodal fermions in a d-wave superconductor are de- 
scribed by {H 

Sf = J d 2 x dr xHd T - A t t z — v F r z id x + vfA x 

-A a T na + VfA^ftffr, - v a t s {e ls * } (id y )] X (Bl) 
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where the electron operators Xaa are defined as 



+ J^drz - A%) + ,r x {z^d x z - Al) J (B6) 



X21 




c- 

k F +k\ 
-kp-ki 


Xl2 
X22 _ 




C k F +kl 
L — fcp — fe| -1 



(B2) 



and the coordinate system was rotated in such a way 
that the x axis goes along the nodal direction that we 
are considering ( see Figure]^). 




FIG. 7. Order parameter for a d-wave superconductor. 
Gapless excitations exist at k F = (±fc_p,0), (0, ±k F ). 

Antifcrromagnetic fluctuations are introduced via 
Saf = J dr J d 2 kd 2 qn ? cl^ + ^ a a al3 c SF+u ^ + h.c. 

d 2 x dr n(x)xaa{x)e ab e^a 7 pXbf3{x) (B3) 



where ft = (cos 0, sin 0). 

Spin and charge may again be decoupled by rotating 
the fermions as in (0) 



(B4) 



Xaa e 

with the result 

5 = y d 2 xdr ijj^[d T - v F T z id x — v^T x {id y )]4> 
+ J d 2 xdri}j aa {x) e ab e^<7*pip b f3(x) 

+ ~J d 2 xdr f^[T z d T (p-2A c T T z - v F d x (p + 2v F A x }^ 
+ ^[a z d T 9 - 2A°t z - v F T z cr z d x e + 2v f t z o z A° x \^ 



This describes the same coupling of the quasiparticle cur- 
rents to the fluctuations of charge and spin as in ( |50| ) 

S f =S° f + J d 2 x dr ( JC(d T (p - Ar) + J°{d x y - A x ) 



with 



J 



- ^a z 4> 



J. i 



-V F tjj^tp 

-v F i\r <j z t z i\i 



(B7) 



Quantum disordering of the superconducting and anti- 



fcrromagnetic orders in (B6) may now be achieved by 



condensing vortices and merons with the possibility of 



five phases similar to phases 3a-3e in Section IV D 



A) Spinons and holons confined. No quantum number 
separation. 

B) Spinons unbound and holons glued to fermions. 

C) Holons free and spinons bound to fermions. 

D) Spinons and holons bound together, decoupled from 
fermions. 

E) All excitations decoupled. Free holons, spinons, neu- 
tral fermions. 



APPENDIX C: ADIABATIC CONTINUATION 
BETWEEN DIFFERENT PHASES OF A Z 2 x Z 2 
GAUGE MODEL WITH MATTER FIELDS. 

In the pure x Z2 gauge theory there are 5 phases 
that are distinct and separated by phase transitions. A 
question that we address in this section is whether this 
distinction survives in the presence of matter fields. 



1. Toy Models 

Let us begin with a simple model 

= -Ki n °y ~ K2 J2 n ■ ~ 13 zZ a a T io v i v i ( ci ) 



where m = ±1 is an Ising matter field. To construct 
the full phase digram of this model we consider several 
limiting cases. When (3 = we have two independent 
gauge fields, each of which has confining and deconfin- 
ing phases. When both K\ and K2 are small we have 
a confining phase for both a and r, that has no extra 
degeneracy on topologically non-trivial manifolds and is 
labeled 1 on Fig. || When K\ is large and K% is small we 
have a phase that is confining for r, but deconfining for a 
Jjihase 2' in Fig. ||). There is an analogous phase for if 2 
v Targe and K\ small (phase 2" in Fig. ||) which is confin- 
ing for a and deconfining for r. When both K s are large, 
we have a fully deconfining phase with degeneracy 4 on a 
cylinder. When K\ — 00 (BCGF plane on fig |J) there are 
no frustrated plaquettes for er, so we can choose a gauge 
where all er,, 



1. The model is then the same as in 



24 



and its phase diagram can be easily constructed. When 
K\ = (ADHE plane on fig ||) we find that integrating 
out <Jij and Vi only adds a constant to the action for r and 
does not affect the confinement-deconfinement transition 
which takes place for the same value of K2, regardless 
of the value of (3. When — 00 (EFGH plane on fig ||) 
we must have Yin <J i3 T n = 1 on every plaquette, so we 
can choose a gauge where &ijTij = 1 on every link. The 
fields a and r are identical and there are only two phases, 
a confining (phase 1) and a deconfming (phase 2) with 
the transition determined by K\ + The full phase 
diagram may now be obtained by connecting the lines on 
the faces of the cube on Fig ||. It is immediately clear 
from this picture that the two partially confining phases 
which appeared to be distinct for (3 — (phases 2' and 
2" in the ABCD plane) may be continuously connected 
through a path that takes advantage of the finite (3 region 
of the phase diagram. It is important to realize that our 
argument for the existence of a path connecting phases 2' 
and 2" does not depend on the details of how the phase 
boundaries in Fig || are connected. One can always find 
a path which begins in phase 2', approaches face ABFE, 
goes up to EFGH, crosses to EHDA, and finally comes 
down to 2" without crossing the phase boundaries (this 
path does not have to actually be on any of the faces 
and it may be sufficient to be in their vicinity). It is in- 
teresting to note that the cross-section DBFH on Fig || 
the phase diagram looks similar to a liquid-gas phase di- 
agram, where the two phases 2' and 2" may be separated 
by a lst-order transition or continuously connected arond 
the critical point which terminates the lst-order line. 
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FIG. 8. Phase diagram of (CI). Plane ABCD corresponds 
to P = 0, EFGH to /3 = 00, ADHE to K x = 0, BCFG to 
K x = 00, ABFE to K 2 = 0, and DCGH to K 2 = 00. The 
pure gauge model (/3 = 0) has only two partially confining 
phases (2' and 2") because not all the possible lattice Maxwell 
terms are present. They are separated by one first order or 
two second order phase transitions for small /3, but may be 
continuously connected for larger /3. 



theory ( pl| ) may be connected to each other is that both 
of them are related to the Higgs phase (phase 2) for the 
Vi matter field. This may be explained by noting that in 
such Higgs phase visons of either a or r are forbidden, but 
their composite is not, so this phase should be related to 
the phases where these visons are condensed separately 
(but not simultaneously). 

Note that this argument may no longer apply if the 
matter field carries a quantum number, and a Higgs 
phase breaks some continuous symmetry. Let us, for ex- 
ample, explore the model where the matter field is an 
XY order parameter 



5 = -^ELK-^£rh: 

□ □ 



□ □ 
i j cos ( 



(C2) 



In this case, the Higgs phase has superfluid order, and 
therefore is fundamental ly distinct from confined insulat- 
ing phases. Thus, we can no longer easily claim the equiv- 
alence of the two phases in which either a or r (though 
not both) fields are confining. 

We can again attempt to construct a phase diagram 
following construction on each of the outside faces of the 
cube. The ABCD plane is the same as in Fig ||. The 
BCGF plane (K\ = 00) will now have three phases: a 
confining and a deconfming phases without broken XY 
symmetry (phases 2" an 4), and a phase with broken 
XY symmetry f§. When K x = (ADEH plane) we 
have 4 phases. This is obvious from the fact that when 
we integrate out <xy we find that Ty and cos4>i are de- 
coupled from each other, and we have separate order- 
disorder and confinement-deconfinement transitions. We 
therefore find two superfluid phases SFi and SF2 that 
differ in their degeneracy on the nontrivial manifolds (on 
the cylinder it is 1 for SFi and 2 for SF 2 ). The origin 
of this extra degeneracy for SF2 is that it has a finite 
energy topological excitation: a bound state of a and t 
visons that does not interact with the matter field. This 
has interesting implication that we have hc/2e vortices 
that are bound to either a or r visons, and the two kinds 
of vortices are distinct. 

We do not at this stage know what the generic phase 
diagram for ( |C2] ) in the K\K%fi cube is. One possibility 
is shown on Figure ^|. As in the Ising case there is a 
way to connect phases 2' and 2" continuously by going 
to finite f3. 



The real reason why phases 2' and 2" of the gauge 
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FIG. 9. One scenario for the phase diagram of (02) with 
two kinds of superfiuid phases SF\ and 3F%. Shaded fig- 
ure shows a phase boundary of the superfiuid and insulating 
phases. Diagonal shading corresponds to the boundary of SF\ 
and horizontal shading to a boundary of SF^.. There is a con- 
tinuous path to go between partially confining phases 2' and 
2" without crossing the phase boundaries. 

There is another qualitatively different phase diagram 
(without introducing new phases) where the point M is 
at the phase boundary with the superfiuid phase. This 
will remove the possibility of a continuous path between 
phases 2' and 2" (see Figure |l0| ). 
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FIG. 10. Another scenario for a phase diagram of (C2) 
when the point M is exactly on the superfluid-insulator phase 
boundary. As before shaded figure shows a phase boundary of 
the superfiuid and insulating phases with diagonal and hori- 
zontal striping that correspond to SFi and SF2 respectively. 
In this case there is no continuous path connecting phases 2' 
and 2" without crossing the phase boundaries. 
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At this point we are unable to make a definite comment 
on of the validity of either of the scenarios shown on 
Figure ^ or Figure [To]. We note however that this issue is 
amenable to study by numerical or other means. Thus, 
future work should be able to settle this satisfactorily. 

Another important model to consider is one in which 
the matter field is fermionic. An appropriate model is 

s = -k y n a v - k ? n r y ~ 13 n a *3 T v^i ( c3 ) 



where the ^'s are real fermions. Following the same kind 
of arguments as before we find the phase diagram shown 
on Fig[n| There is no Higgs phase for the fermions which 
leads to phases 2' and 2" being distinct even for finite (3. 



FIG. 11. Phase diagram of (C3). When the matter field is 



fermionic distinction between phases 2' and 2" survives for all 
/3. 



2. Full Action 

Let us now consider the action ( |105| ) and ask how 
many truly distinct phases it has. The phase space of 
this model is large and an explicit construction of the 
full phase diagram is difficult. We note however that 
the charge sector o f the t heor y is precisely the same as 
(|C2|). Consequently if in |C2j , the two phases 2 and 2' 
are smoothly connected to each other, they will neces- 
sarily be so f or the full action as well. If on the other 
hand, in (C2), the two phases are distinct, then that is 
evidence (though not proof) that they are distinct in the 
full theory as well. 

Thus, unambiguous determination of the phase dia- 
gram of (C2) will shed considerable light on the impor- 
tant conceptual issue of whether CbSf and CfSb are dis- 
tinct or not. 

While the distinction between CbSf and CfSb, if any, 
is subtle, it is very clear that they are both distinct from 
CSbNf. One can not find any vison attachment that 
would map the spectrum of CSbNf to either CbSf or 
C fSb, which proves rigorously that it is a phase funda- 
mentally distinct from the other two. 

To summarize the discussion in this section, we argued 
that studies, nume rica l or otherwise, of simple models of 
the form of Eqn. (C2) should be extremely useful in de- 
ciding on the issue of whether CbSf and CfSb are distinct 
quantum phases. One can, however, prove rigorously 
that CbSf and CfSb are fundmenta lly d ifferent from the 
other partially c onfin ing phase of (105) CSbNf. The 
other phases of (105): the fully confining phase (CSf) 
and the fully deconfming phase (CbSbNf) will be dis- 
tinct from any of the partially confining ones and from 
each other as may be seen from their degeneracy on non- 
trivial manifolds. 



APPENDIX D: UN-FRACTION ALIZED PHASES 

In this paper, we have, for the most part, focussed 
on states in which the electron is fractionalized. How- 



2G 



ever, even the transitions which do not lead to electron 
fractionalization are rather novel. One would ordinar- 
ily assume that strong quantum fluctuations will com- 
pletely disorder a p-wave superconductor. However, as 



we pointed out in section IV B, if hc/Ae vortex- 7r disclina- 
tion composites are gapped, then the spin symmetry can 
be restored without affecting the charge; alternatively, 
the superconductivity can be destroyed without affecting 
the spin ordering. 

Let us consider, first, what happens when flux hc/2e 
vortices condense, but no other topological defects con- 
dense. Then the charged degrees of freedom are disor- 
dered, but the spin nematic order parameter should be 
undisturbed. Following the arguments of pij , a possible 
un-fractionalized spin nematic insulating state is (in the 
notation of ff3|| ) a triplet p x density wave: 

(</> Qt (fc + Q,i)<MM)) =®Q-d% sink x a (Dl) 

This state is related to the p-wave superconducting state 
by a 'rotation' generated by 



d 2 k 



z\(k) c\(-k + Q) 



(D2) 



In other words, the triplet p x density wave and the p- 
wave superconductor are related in precisely the same 
way as a CDW and an s-wave superconductor. In par- 
ticular, Hamiltonians with short-ranged interactions can 
be constructed for which both states are exactly degen- 
erate; such Hamiltonians could describe a critical point 
between these two states. 

In the triplet p x density wave state, there is no spin 
moment (at any wavevector), since the right-hand-side 



of (DI) vanishes upon integration over k. However, the 
spin-nematic order parameter, which may be calculated 
from (Dl), is non- vanishing: 



l s » s2 



diag(2/3,-l/3.-l/3) 

(D3) 



where Si is the i-th component of the total spin of the 
system. Hence, this is the natural spin nematic state 
which results when a p-wave superconductor is quantum 
disordered by flux-/ic/2e vortex condensation. The pos- 
sibility of spin nematic states in the context of high Tc 
cuprates has been proposed in |47j . 

When the spin degrees of freedom are disordered, but 
the charge remains ordered, the triplet p-wave super- 
conducting order parameter and the spin-nematic or- 
der parameter vanishes; only the charge-4e order pa- 
rameter is left. The condensation of merons causes the 
fermionic quasiparticles to be confined to spin. Once the 
merons have condensed, the topological quantum num- 
ber in the spin sector is no longer well-defined, so the 



flux hc/Ae vortex- tt disclination composites become sim- 
ple flux hc/Ae vortices, as we would expect for a charge- 
4e superconductor. Said differently, the meron conden- 
sate screens the spin topological charge of the flux hc/Ae 
vortex- 7r disclination composites, thereby making render- 
ing them simple flux hc/Ae vortices. Remarkably, by 
quantum disordering the spin sector of a p-wave super- 
conductor, we have changed the charge of its order pa- 
rameter, which may, for example, lead to some unusual 
critical behavior of the superconducting transition. 
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